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ABSTRACT 

We analyze the global structure of 34 late-type, edge-on, undisturbed, disk galaxies spanning a wide 
range of mass. We measure structural parameters for the galaxies using two-dimensional least-squares 
fitting to our i?-band photometry. The fits require both a thick and a thin disk to adequately fit 
the data. The thick disks have larger scale heights and longer scale lengths than the embedded thin 
disks, by factors of ~2 and ^^1.25, respectively. The observed structural parameters agree well with 
the properties of thick and thin disks derived from star counts in the Milky Way and from resolved 
stellar populations in nearby galaxies. We find that massive galaxies' luminosities are dominated by 
the thin disk. However, in low mass galaxies (Vc < 120 km s~^), thick disk stars contribute nearly 
half of the luminosity and dominate the stellar mass. Thus, although low mass dwarf galaxies appear 
blue, the majority of their stars are probably quite old. 

Our data arc most easily explained by a formation scenario where the thick disk is assembled through 
direct accretion of stellar material from merging satellites while the thin disk is formed from accreted 
gas. The baryonic fraction in the thin disk therefore constrains the gas-richness of the merging pre- 
galactic fragments. If we include the mass in HI as part of the thin disk, the thick disk contains < 10% 
of the baryons in high mass galaxies, and ~ 25 — 30% of the baryons in low-mass galaxies. Our data 
therefore indicate that the fragments were quite gas rich at the time of merging {fgas = 75 — 90%). 
However, because low mass galaxies have a smaller fraction of baryons in their thin disks, the pre- 
galactic fragments from which they assembled must have been systematically more gas poor. We 
believe this trend results from increased outflow due to supernova-driven winds in the lower mass 
pre-galactic fragments. We estimate that ^ 60% of the total baryonic mass in these systems was lost 
due to outflows. Pushing the episode of significant winds to early times allows the mass-metallicity 
relationship for disks to be established early, before the main disk is assembled, and obviates the 
difficulty in driving winds from diffuse disks with low star formation efficiencies. We discuss other 
implications of this scenario for solving the G-dwarf problem, for predicting abundance trends in 
thick disks, and for removing discrepancies between semi-analytic galaxy formation models and the 
observed colors of low mass galaxies. 

Subject headings: galaxies: formation — galaxies: structure 



1. INTRODUCTION 

The structure of galactic disks provides strong con- 
straints on their formation and evolution. Spiral galax- 
ies have long been recognized to contain several distinct 
populations of stars (e.g., disks, bulges and halos), each 
with distinct chemical and kinematic properties that cap- 
ture unique epochs in the formation of the galaxy. Ob- 
servations of the Milky Way and a wide range of other 
galaxies have revealed the need for yet another compo- 
nent, namely, a thick stellar disk. Originally detected as 
an excess of light at high galactic latitudes in deep sur- 
face photometr y of early-type galaxies (Burstcin 197^ 
iTsikoudi 1979), a thick disk was later revealed in the 
Milky Way using star counts l)Gilmore fc Reid 1983') . 

The properties of the Milky Way's thick disk have 
revealed many differences from the thin disk. Struc- 
turally, the Milky Way's thick disk has a significantly 
larger scale height t han the thin disk, as its name im- 
plies ( for reviews se e lReid fc Ma'iewski. 1993; .Buser et al.l 
119991 iNorrisI 119991 and references therei n) . It also 
may have a somewhat longer scale length (jRobin et al.l 

^ e-mail address: 

^ Alfred P. Sloan Foundation Fellow 
^ e-mail address: 



T99i lOihall200 'it lChen et a,ljl2 00Tt ILarsen fc HnmnhrevFl 



2003j) . Thick disk stars are older and more metal- 
poor than stars in the thin disk (e.g. iReid fc Maiewskil 
119931 Ghiba & Beers 2000). They are also signif- 
icantly enhanced in a-elements, compare d to thin 
disk s tars of comparable ir o n abundance l | Fuhi 
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Bensbv et al. "2003"; 'Feltzine et al. "20G3!;'Mishenina et ajj 
1)04; Brewer & Carnev 2004; Bensbv et al. 2005). Kine- 
matically, Milky Way thick disk stars have both larger 
velocity dispersions and slower net rotation than stars 
in the thin disk pissml IT995t iChiba fc Beer j l200?i 
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For many years, however, it remained unclear whether 
the thick disk was a truly distinct component of the 
Milky Way, or whether it was only an older, metal- 
poor extension of the thin disk, as might be created by 
stead y vertical heating over t he lifetime of the Galaxy 
(e.g., iDove fc Thro nson 1993). Over the past five years, 
conclusive evidence that the thick disk is indeed distinct 
from the thin disk has come from a series of detailed 
chemical abundance studies. Stars with thick disk kine- 
matics show significant alpha-enhancement compared to 
thin disk stars with identical iron abundances, thus form- 
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ing a separate parallel sequence i n a plot o f fa/H] vs 
[Fe/H] (see the recent review bv iFeltzing et al. 20041). 
Studies of resolved stars in nearby galaxies also find a 
thick disk of old red giant branch stars whose lack of 
metallicity gradient cannot be explained by steady ver- 
tical heating ( Seth et al. 2005. : Mould 200i . 

Three general classes of formation mechanisms have 
been proposed to explain the properties of the Milky Way 
thick disk. In the first, a previously thin disk is dynam- 
ically heated to form a thi ck disk, afte r which a new 
thin disk forms llQuinn et al. 199a; Velazauez k White! 
^99; Robin et al 199fit iChen et al.ll2001() . In the sec- 
ond, the thick disk forms directly from gas at a large 
scale height, possibl y during a largely monolithic proto- 
galactic co llapse (Eggen et al. 1962; Gilmore & Wvse 
|1986t INorris & Ryan 1 991; Burkert et al. 1992; Kroupa 
12002; Fuhrmann H20b4l: iBrook et al.li200ij) . In the third, 
the thick disk forms from a series of minor-merger 
events which directly deposit sta rs at large scale heights 
l|StatleillT988l lAbadi et alJl2003j) . Recent cosmological 
simulations have suggested a more comp l icated origin. 
Disk galaxy simulations by lAbadi et al.l l)2003j) find a 
thick disk which is composed primarily of tidal debris 
fro m disrupted s atellit es while comparable simulations 
by IBrook et al.l l)2004|) find that thick disk stars form 
during a period of chaotic mergers of gas-rich building 
blocks. Recent kinematic measurements favor scenarios 
wher e mergers play a significant role in thick disk form a- 
tion l|Gilmore et alJl2002t lYoachim fc Dalcantonll200,5t) . 

While all of the above scenarios are viable expla- 
nations for the origin of the Milky Way, the struc- 
tural parameters of thin and thick disk components in 
a wide range of galaxies can help distinguish among 
these formation scenarios. Unfortunately, the mea- 
surements required to characterize thick disks are diffi- 
cult to make outside the Milky Way. The Milky Way 
thick disk provides less than 10% of the local stel- 
lar density ijBuser et al.lll99^ . and this faintness ham- 
pers detailed study of comparable extragalactic thick 
disks. To date, thick disk structural properties have 
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de Griis & Peletieill997l see Table 2 below). These stud- 



ies analyze galaxies in the edge-on orientation, which al- 
lows clear delineation between regions where thin and 
thick disk stars dominate the flux. The edge-on orienta- 
tion also allows line of sight integrations of faint stellar 
populations to reach detectable levels. 



In this paper, we analyze a large sample of edge-on 
galaxies and decompose them into thick and thin disk 
components. Analysis of i?, R, & Ks photometry and 
color maps has previously revealed that these galaxies 
are surrounded by a flattened faint red envelope, with 
properties very similar to th e Milky Way thick disk 
(jDalcanton & Bernstein 2002^) . We now use a full 2- 
dimcnsional fitting procedure capable of simultaneously 
fitting the thick and thin disk light distributions to derive 
their full structural parameters. 

1.1. Galaxy Sample 

The sample used in this paper was drawn from the opti- 
cal and infrared imaging found in lDalcanton fc Bernstei'nl 



TABLE 1 

Galaxies rejected from the 2-disk H-band fitting 



FGC(E) 


Reason for rejection 


51 


Udliie ^^eilLl dl JvliOL 


84 


kJIJIidl d± lilt!) vlOIUlt. 


143 


Bright central star-forming region 


442 


Spiral arms visible 


256 


Thin disk below seeing limit, fits did not converge 


1863 


Large warp, bright foreground stars 


1945 


Fits did not converge, possibly spiral arms 


1971 


Polar ring galaxy 


2217 


Large bulge component 


2367 


Spiral arms visible 


2264 


Fits did not converge, scattered light problem 


2292 


Thin disk below seeing limit, bright foreground stars 


E1440 


Asymmetric disk 


E1447 


No velocity data 


E1619 


Large bulge component 



("2000", hereafter Paper I). Briefly, our sample of edge- 
on bulgeless galaxies was i nitially selected from the F lat 
Galaxy Catalog (FGC) of IKarachentsev et all ljl993|) . a 
catalog of 4455 edge-on galaxies with axial ratios greater 
than 7 and major axis lengths grater than 0.6'. The 
color maps and initial det ections of the thick disk s in 47 
galaxies were presented in lDalcanton fc BernsteinI l)2002l 
hereafter Paper II). 

Not all galaxies from Paper I have been included in 
the analysis presented here. We have excluded several of 
the more massive galaxies with sizable bulge components 
that could not be adequately masked or modeled. We 
likewise eliminated several low mass galaxies with bright 
central star clusters for similar reasons. We have also re- 
moved any galaxies that have either significant warps or 
visible spiral arms (i.e., that were not viewed perfectly 
edge-on), as these systems are poorly modeled by our 
fitting procedure. Finally, we eliminated galaxies whose 
surface brightness profiles would be severely affected by 
atmospheric seeing. A full list of the 15 excluded galax- 
ies are listed in Table leaving a sample of 34 galax- 
ies suitable for decomposing into thick and thin com- 
po nents. When possib l e, we have used distances listed 
in IKarachentsev et al.l lj200Cltl derived from a local flow 
model. Otherwise we use the galaxy's redshi ft correcte d 
for the motion of the Local Group (|Yahil et alJll977j) . 
assuming a Hubble Constant of Hq = 70 km s""'^ Mpc""'^. 

2. 2D FITTING 
2.1. Galaxy Models 

The distinctive vertical color gradients identified in Pa- 
per II suggest that the stellar population above the galax- 
ies' midplanes is different from the one within it. We as- 
sume this change is due to the existence of two distinct 
stellar populations analogous to the MW's thick and thin 
disks. Our 2-dimensional fitting procedure attempts to 
decouple these two populations to measure their scale 
heights, scale lengths, and luminosities. 

We model the surface brightness of each disk compo- 
nent as a radially exponential disk. We adopt the lumi- 
nosity density £ of each disk component to be 

/:(i?,z)=£oe-«/'^«/(z) (1) 
where (i?, z) are cylindrical coordinates, Co is the central 
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luminosity density, hn is the radial scale length, and f{z) 
is a function describing the vertical distribution of stars. 

Throughout, we adopt a generalized vertical distribu- 
tion 



/(z) = sech2/^(7Vz/zo) 



(2) 



where zq is the vertical scale height and A'^ is a parameter 
controlling the shape of the profile near the midplane. 
For appropriate choices of N, this equation can repro- 
duce many popular choices for the vertical distribution of 
star light. With iV = 1 EquationlUbecomes the expected 
form for a self-gravitating isothe r mal sh eet l)SDitzerll942l: 
lva,n der Kruit Sea,r1elll98m ITflSl. When ^ oo, 
Equation 13 reduces to f{z) cx e"^'"'', where — zo/2. 
Previous fits to the vertical distribution suggest that an 
inter mediate value of = 2 is a better model of galaxy 
disks l)van der Kruitl[T988|l . as expected for the superpo- 
sition of several population s with a range of vert ical ve- 
locity dispersions (dc Grijs fc van der Kruitll99ff> . How- 
ever, different values of N only produce differences near 
the galaxy midplane, and all share exponentially declin- 
ing profiles at large radii. 

When fitting a thick plus thin disk model, we prefer- 
entially use TV = 1 for both components because of its 
physical motivation. We note that since our main goal 
is not to model galaxies near their midplane where these 
functions have their largest differences, our results are 
not particularly sensitive to the choice of model. To per- 
mit comparisons to previous work, we also derive single 
disk fits to our sample using Equation |21 without a fixed 
N, allowing the shape of the vertical profile to vary to 
best fit the data. 

To translate the adopted luminosity density into the 
observed surface brightness distribution, we assume that 
the disks are viewed perfectly edge-on. Other au- 
thors have demonstrated that slight deviations from 
i — 90° have minimal impact on the derived struc- 
tural parameters (e .g., Ivan der Kruit fc Searlel Il981bt 
Ide Griis et al.l Il997t) . We also assume that scale 
heights are independent of projected radius for late-type 
galaxies , as found by 'van der Kruit & Scarlc (1981b); 



Bizvae v & Mitronova (2002): Ide Griis fc P clcticr (1997); 
Shaw fc Gilmore (^1990). With the above assumptions. 



(3) 



the model edge-on disk surface brightness is given by 

I](i?,z) = So,o(i?//ii?)A'i(i?/M/(^) 
where K\ is a modified Bessel function of the first or- 
der, I]o,o is the edge-on peak surface brightness (So,o = 
2/ifl;£o), and R is now the projected radius along the 
major axis. The face-on surface brightness of such a 
disk is E(i?) = Soe"^/'* with Eo = 2zo/:o- Through- 
out, we convert our edge-on peak surface brightnesses 
to magnitudes using ^(0, 0) = rtizp — 2.51og(Eo,o) where 
vfizp is the photometric zero point from Paper I. The 
face-on central surface brightness can then be calculated 
as /xo = ^^(0,0) — 2.51og(zo//ifl)- The conversion be- 
tween the edge-on and face-on orientation assumes that 
disks are optically thin at any orientation, an assump- 
tion that is obviously not true for massive galaxies with 
dust lanes. However, we correct for this effect later in 
ii4.3l We do not model any possible disk truncation, as 
this is a small effect seen only in the region R > S/tp 
(Ivan der Kruit fc Searle Il981bl: iKreeel fc van der KruitI 
l2004HPohlen et alJl200flt) . 



Our sample of galaxies was initially selected to be 
"pure disk" systems, and thus there are very few galaxies 
which possess a prominent bulge component. We there- 
fore do not attempt to decompose a bulge component 
from the surface brightness distribution and simply re- 
ject galaxies with significant bulges from the sample (Ta- 
ble [fl). 

We have tested if the profiles described by Equation [5] 
could be significantly affected by seeing. We convolved 
model images with a two-dimensional circular Gaussian 
kernel to simulate the atmosphere's effect. We found 
that this step, in general, was unnecessary. Unconvolved 
fits differed from convolved fits only for the most distant 
galaxies. Several of these galaxies have been eliminated 
from the sample, as listed in Tabled 

2.2. Fitting Method 

We use Levenberg-Marquardt least squares fitting of 
the galaxy images to find the best parameters for the 
models described in H2.ll Before fitting, the images of 
the galaxies are sky-subtracted and foreground stars and 
background galaxies are generously masked (see Paper 
H). The images are cropped at i? ~ Aha to speed com- 
putation time. Our tests have shown that the fits are 
insensitive to the exact cropping, causing variations in 
individual parameters of only a few percent. The crop- 
ping also reduces the chance that our fits could be biased 
by warps or flaring of the disk s at large radi i. 

Following the technique of IKreeel et alJ ()2002|) . we 
weight each pixel by the inverse of the model surface 
brightness distribution at that pixel. By using the model 
rather than the data to determine the weighting, we elim- 
inate the bias of overweighting positive noise spikes. This 
weighting scheme places large amounts of weight on the 
lowest signal-to-noise pixels, ensuring that regions of low 
surface brightness (i.e., where a faint thick disk could be 
detected) receive adequate weighting. To prevent the fit 
from being overwhelmed by regions with low signal-to- 
noise, we set the weight to zero beyond the 1-a noise 
contour, defined as where the model falls below the stan- 
dard deviation of the background. Due to the low signal- 
to-noise ratio in the band data, these images were 
clipped at the 1/2-a level to ensure that an adequate 
number of pixels were included in the fit. Each fit was 
iterated four times to ensure convergence of the model 
parameters and weighting scheme. Fits were performed 
using pixel coordinates and counts, then converted to 
arcseconds and magnitudes using the calibrations in Pa- 
per I. 

It is common practice when fitting models to edge- 
on galaxies to cro p out regions nea r the midplane 
of the disk (e.e., Kregelet alJ l20f)l Ide .Tond ll9M 
iBizvaev fc Mitronova.2002,) . Cropping avoids the hard- 
to-model effects of dust lanes, bulges, and star forming 
regions. The color maps of our galaxies imply that galax- 
ies rotating at speeds less than 12 km s~^ do not c ontain 
concentrated central dust lanes (jDalcanton et al.l l2004'l. 
For more massive galaxies, our weighting procedure en- 
sures that any midplane structure receives a minimal 
amount of weighting when calculating the goodness-of-fit 
X^- We chose to fit models both with and without the 
midplane cropping to quantify the systematic uncertain- 
ties introduced by midplane structure. 

We begin by fitting single disk models to all three B, 
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R, and, Ks images, holding the galaxy position and ro- 
tation fixed. We then fit 2-component models to the 
images, allowing the ofi^sets and rotation to vary, but 
constraining all components to have the same center and 
orientation. For this second step, we use only the R- 
band images due to their high signal-to-noise. Ideally, 
we would perform the 2-disk decomposition in the Ks 
band which best represents the smooth stellar distribu- 
tion and is least effected by dust. However, due to the 
bright infrared sky, the NIR images are lower signal-to- 
noise and cannot reach to faint regions where a thick 
disk would dominate. The i?-band therefore represents 
the best compromise between reaching faint regions of 
the galaxies while minimizing the effects of dust extinc- 
tion and bright star forming regions. 

When fitting a single disk with only three free param- 
eters, our procedure converges to the same minima 
given any reasonable initial guesses. However, for multi- 
ple component models, which have up to 10 free param- 
eters, we find that fits often converge to local minima 
rather than to the global minimum. To ensure we find 
the global minimum when fitting multiple components, 
we fit each galaxy using up to 50 un ique initial parameter 
guesses, following IWu et all pOO^ . The initial parame- 
ters for each galaxy model were randomly varied up to 
±50% to ensure we cast a large net in parameter space. 

The formal parameter uncertainties that result from 
our fits are not meaningful because we used a weighting 
scheme that is not based on the actual pixel uncertain- 
ties. Even if we did minimize using formal pixel errors, 
our returned uncertainties would be much too low. The 
formalism requires residuals to be Gaussian, which is 
rarely the case when fitting nearby galaxies. The situ- 
ation is comparable to trying to model Mount Rainier 
as a cone-you can do it, but the residuals will be domi- 
nated by real physical structures and not random Gaus- 
sian measurement errors. In the case of spiral galaxies, 
real substructure exists in the form of spiral arms, dust 
lanes, regions of active star formation, warps, fiares, HII 
regions, etc. As an alternative assessment of the system- 
atic errors which are likely to dominate our uncertain- 
ties, we fit a series of models using a variety of different 
weighting and masking schemes (Table 4) and quote the 
median result for each parameter. We then adopt the 
full range of convergent models for each parameter as a 
measure of the inherent systematic uncertainties. The 
resulting uncertainties are 2-100 times greater than our 
formal x^ uncertainties, confirming that systematic er- 
rors dominate our uncertainties. 

2.3. Tests on Artificial Images 

To assess the reliability of the 2-disk decompositions, 
we created a set of 100 artificial galaxies. We adopted the 
surface brightness profile in Equation |2 with an = 1 
vertical distribution for both a thick and thin component, 
and varied the structural parameters of the disks within 
ranges similar to our sample galaxies (for the thin disk: 
20.7 mag/n"< ^(0,0) < 22.7 mag/D", 2.5"< Hr < 
19.4", 0.6"< zq < 3.6"; and the thick disk had parame- 
ters in the range: 21.4 mag/n"< fj.{0, 0) < 24.2 mag/D", 
1.5"< hR < 51.2", 1.4"< zo < 20"). The model galaxies 
were convolved with a circular Gaussian with 1" FWHM, 
typical of the seeing for the observations. Wc then added 
read noise, sky noise, and Poisson noise to the simulated 



galaxies, with amplitudes chosen to mimic our _R-band 
data. The galaxy images were rotated up to 2 degrees 
and offset up to 2 pixels (0.5") from the image center. We 
then fit the galaxies with the two disk models, with and 
without seeing corrections. We assume = 1 and use 
the same spread of initial parameter guesses as described 
in 

Of the 100 simulated galaxies, only three fits failed to 
converge. 90% of the scale lengths are recovered to within 
±2% of the input value, with all of the results converg- 
ing within ±10%. 90% of the scale heights are recovered 
within ±3% of the input value with all results within 
±12%. 90% of the central surface brightnesses converge 
to within 0.09 mags of the correct result. The orienta- 
tions were always correct to within 0.1 degrees, with a 
median error less than 1%. All of the spatial offsets were 
within 1 pixel of the correct position. There were no 
systematic trends in the size of the errors verses galaxy 
properties. The high accuracy of these fits indicates that 
we are not limited by pixel noise in our fits. However, 
since our model was a perfect match to the input data, 
this correspondence is not surprising. 

We tested models that did not make a seeing correction 
convolution and found the fits still returned scale lengths 
and heights that were accurate to < 0.2", as long as zq > 
1". The majority of our observed galaxies do have zq > 
1", and thus we do not account for seeing i n our fits. This 
result is consistent with the analysis of Ide Griis et alJ 
(1997) who find that for an exponential vertical profile, 
convolution is unnecessary when the seeing FWHM < 
0.6h,. 

In addition to testing our ability to recover the pa- 
rameters of a known model, we also tested our ability to 
correctly measure the structural parameters when using 
an incorrect function for the vertical light distribution. 
Specifically, we fit an A^ = 2 model instead of the correct 
A^ = 1 vertical profile to each disk. These fits returned 
results similar to the fits using the correct model, and the 
resulting scale heights and lengths fell within ~ 5% of the 
correct values. This indicates that the galaxy sizes are 
constrained primarily by light well away from the mid- 
plane. The luminosities were more divergent, however, 
due to the large differences between these models at their 
midplanes. The N ~ 2 model was slightly biased towards 
having over-luminous thick disks (with a few outliers as 
well), but the majority (70%) of fits were within a factor 
of 2 of the correct Lthick/ Lthim despite being fit with the 
wrong function. 

Finally, we also tested our two-disk fitting code on arti- 
ficial galaxies that had no second thick disk component. 
In these cases, the fits always converged to extremely 
faint thick disks (< 1% of the thin disk flux) and usu- 
ally converged to either very large or very small thick disk 
scale lengths, mimicking either a uniform sky background 
or a small point source. Overall, these results encourage 
us to believe that if there are no thick disk components in 
our data our fitted parameters will diverge to unphysical 
values. 

3. SINGLE DISK FITS 

Before discussing the results of decomposing the galax- 
ies into two disk components (fQ, we discuss the results 
for fitting single disks to the light distributions. These 
fits are useful simple descriptions of the galaxies, and the 
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galaxy circular velocity. In Figure 13 we compare our sin- 
gle disk R-an d -B-band fitted scale lengths with the edge- 
on sample of 
surements of 



Kreeel et all EOO: 



Fig. 1. — Cumulative distributions showing the systematic effects 
of fitting our galaxies with varying models. All models are com- 
pared to a fit using inverse weighting and no masked regions. Fits 
for all three filters have been combined. The solid curve shows the 
model which used uniform pixel weighting, while the dotted curve 
shows the midplane masked model, and the dashed curve shows 
the \R\ < 3/iij radially masked model. 



resulting parameters can be directly compared to previ- 
ous fits of edge-on and face-on galaxies. 

As discussed in > i2.'2\ we quantify our systematic uncer- 
tainties using Equation with an = 1 vertical profile 
and a variety of weighting and masking techniques result- 
ing in five different fits for each galaxy image. The five 
fits are: (1) the full galaxy with inverse model weight- 
ing as described in H2.2I (2) the full galaxy with uni- 
form weighting to more heavily weight the high signal- 
to-noise regions; (3) inverse weighting with the midplane 
region (z = ±zo) masked; (4) inverse weighting with 
the outer region R > 3/i/j masked, to eliminate regions 
where our fit may be affected by unmodeled stellar trun- 
cation; and (5) inverse weighting with the high latitude 
region z > 2zo masked to minimize the effect of thick 
disks. The results of the fits are given in Table 3 for all 
3 band passes. The columns show the median edge-on 
peak surface brightnesses (/i(0,0)), radial scale lengths 
{hn), and vertical scale parameters (zq) along with their 
uncertainties. We emphasize again that these are not 
the formal statistical uncertainties (which are deceptively 
small) , but instead are the full range of values to which 
the five different fits converged. 

In addition to quantifying our uncertainties, the five 
different fits provide insight into how variations in fitting 
methods and weighting schemes affect our results. The 
systematic effects of the different methods are plotted in 
Figure n The most notable features in Figure ^ are the 
large systematic shifts in the values of zq for the single 
disk fits. Models that are weighted to fit the midplane 
(e.g. the uniform weighting model), return thinner disks 
while models that mask the midplane return larger val- 
ues of Zq. This effect is present in all three filters, and 
is exactly what one would expect if disk galaxies were 
dominated by thin disks at their midplanes and by thick 
disks with larger scale-heights in the fainter regions. The 
radially cropped and midplane cropped models result in 
fits that have fainter central regions and slightly larger 
scale heights. This is a strong indication that most of our 
galaxies do not have dust lanes which need to be masked. 
Cropping regions at high z has a minimal (< 5% change) 
effect on the fit parameters. 

The parameters for the single disk fits listed in Table 3 
are plotted in Figures EJ El El El and El as a function of 
the galaxies' circular velocity. In Figure |21 we see the 
expected trends that more massive galaxies have larger 
scale heig hts. We also p lot scale heights from the edge-on 
sample of lKreeel et al.l ( ^2002 ) and find that both studies 
give consistent results for the scale height as a function of 



and the fac e-on mea- 

JMacArthur et aO l|200^ : Ide Jon g (1996), 

and iSwaters fc Balcellsl (120021) . Overall, we find that 
lower mass galaxies in the range 50 < K; < 120 have 
scale lengths wholly consistent with measurements made 
in comparable face-on systems. However, the highest 
mass galaxies in our sample have scale lengths that are 
slightly larger than the average found in previous studies, 
although they are still within the full range of the com- 
parison data. This offset is worse in the B-band than 
in the i?-band, and almost certainly reflects the pres- 
ence of strong dust lanes in the more massive systems. 
The higher attenuation towards the central regions of 
the galaxies will suppress the surface brightness at small 
radii, leading to apparently larger scale lengths. This off- 
set may also explain why studies of edge-on disks suggest 
that disks truncate at only 3-4 hn, whereas face-on stud- 
ies see no obvio us signature of truncation at these radii 
(jBarton fc Thom pson 1997; Wciner et_al. 2001,, but see 
also Pohlen et al. 2002). 

In Figure 21 we compare the structural parameters de- 
rived in different band passes. We confirm that redder fil- 
ters converge to shorter scale lengths (Figure EJ, a result 
of the strong radial color gradients seen in both our sam- 
ple and in face-on gala xies (e.g., iBell & de JondEoOU: 
iMacArthur et al.ll2004|) . 

We also find that for galaxies without dust lanes the 
B-band scale heights are predominantly thinner than 
the i?-band. This offset is consistent with the detection 
of strong vertical color gradients in Paper II, where we 
found that the midplanes of late-type galaxies were typ- 
ically bluer than the light above the plane. Somewhat 
unexpectedly, our i^-band scale heights are also signif- 
icantly thinner than the i?-band. We believe that this 
is due to three effects. First, the if-band data does not 
reach as deep as the other filters, making it insensitive 
to the extended thick component. Second, the thinner 
K-hand scale height may indicate the presence of dust 
which blocks light from the midplane in optical filters. 
Finally, there is some indication from studies of resolved 
stars in nearby galaxies that the K-hand light is not com- 
pletely dominated by old red giant stars, but instead has 
a significant contribution from young stars with small 
scale heights (Scth ct al. 2005). 

In Figure El we plot the axial ratios of our sample 
galaxies. O yerall, our axial ratio s are consistent with 
the work of iBizvaev fc Mitronoval (|^002) who measure 
flatness parameters for 153 edge-on galaxies imaged in 
the 2MASS survey and find values of hr/zo ranging from 
~ 2 to '--^ 10. O ur results are also consistent with the 
axial ratios from iKregel et al.l p0021 . We see a slight 
trend for more massive galaxies to be flatter than less 
massive galaxies. Other studies have also suggested 
that low mass, low surface brightness dwarf galaxies are 
thicker than regular spirals. Estimates of the intrin- 
sic axial ratios of dwarf irr egulars range fr om b/a ^ 
0.3 (jHodge fc Hitchco'cSll966 : Binggcli & Pop escul 19951) 
to b/ a ^ 0.6 ()Sung et all Il998i: iStavelev-Smith et alJ 
1 ^^2j). all of which are rounder than typical spirals (e.g. 
Kudrva et al. "1994V We have discussed possible expla- 
nations for this behavior in Dalcanton et al. (_^2004). 
Figure El shows the edge-on peak surface brightnesses 
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Fig. 2. — Single disk scale heights for the i?-band fits. Galaxies 
with prominent dust lanes are plotted with with open circles. For 
comparis on, we sh ow the /-band scale heights from the edge-on 
sample of lKregeT et al. (2002), plotted as red diamonds. 
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Fig. 3. — Single disk scale lengths for the R-band (left) and 
B-band (right) fits. Galaxies with prominent dust lanes are plot- 
ted with with open circles. For comparison, we show other single 
disk fits gathered from the literature. The Ide Jong ( 1996) and 
IMacArthur et al.l 120031) data are face-on or moderately inclined 
galaxy samples with the scale lengths measured in the R and B- 
ban ds ( plotted as blue squares and green triangles respectively). 
The ISw atcrs & Balcells 1 2002) sample consists of late-type spiral 
and irregular galaxies with scale lengths measured in the R and 
B-band (plotted as red diamonds). The Kregel ct al. 12002) data 
were measured from edge-on galaxies in the /-band (plotted as blue 
squares). 



for the one-disk fits. The peak surface brightness of the 
B-hand data is roughly constant, showing httle trend 
with galaxy mass. However, because the FGC sample 
was initially selected from the POSS-II survey plates, we 
would not expect the S-band surface brightnesses to be 
below /X ~ 23 mag/D". On the brighter end, the Free- 
man law ijFreemaiJ Il970l|) suggests a maximum surface 
brightness for edge-on disks. Thus, the B-band peak sur- 
face brightnesses must be confined to a limited range. In 
contrast, we do see increasingly strong trends of surface 
brightness with mass in the redder filters, and particu- 
larly in Kg- Because the selection criteria for the FGC 
limited the range of B surface brightness the observed 
trends in R and K are due to variations in galaxy color 
with mass. As we will discuss in H4.3I extinction from 
dust prevents us from being able to reliably convert the 
edge-on brightnesses to comparable face-on values. 

4. TWO DISK FITS 

4.1. Need for a Second Component 

The traditional signature of thick disks is the presence 
of excess light at high latitudes after subtracting a single 
disk component. To demonstrate the expected excess, 
we subtract the single disk models from the data and 
sum the residuals (inside the 1-cr noise contour) along 
the major axis. The resulting residuals are plotted in 
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Fig. 4. — Comparison of scale lengths and heights for the single 
disk fits in diff'erent bands. Open symbols are used for galaxies 
with prominent dust lanes. 





Fig. 5. — Single disk fits showing the flatness (hn/zo) for each 
band. Open symbols are used for galaxies with prominent dust 
lanes. Dotted lines show the average flatness for a sample of 34 
galaxies in /-band presented in^regel et al. (2002). Dashed lines 
show the average flatness measured measured from a sample of 153 
galaxies from the Revised Flat Galaxy Cata log imaged by 2MASS 
in the /f-band and presented in lBizvaev &: Mitronovai li2002i') 
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Fig. 6. — Edge-on peak surface brightnesses for the single disk 
fits. Open symbols are used for galaxies with prominent dust lanes. 
Points have not been corrected for internal extinction. 



Figured and demonstrate that the single disk fits from 
f|31 systematically leave excess flux at high latitudes for 
all masses of galaxies. We also average the vertical pro- 
file residuals across different galaxy mass ranges and find 
the two disk model is superior to the single disk model 
in all cases. For a single disk model, we can slightly 
improve the fit at high z by allowing the index N to 
vary. However, on average, the absolute value of the 
two-disk model residuals are smaller than the variable 
N model at every height. By collapsing along the ra- 
dial direction, we are assuming that any disk compo- 
nents have a nearly constant scale height with radius, 
as has been found for late type disks in many stud- 
ies llyan der Kruit fc Sear lc 1981b; Bizvacv & Mitronov^ 
'2nn2":' de Criis fc Peletieijll 99711 . 

These tests show that (1) our galaxies are poorly fit 
at large z by the simple sech^ model; (2) by leaving the 
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Fig. 7. — Left: Residuals collapsed along the radial direction 
and averaged over galaxies binned by mass. The single disk models 
(dashed lines) are a very poor fits, leaving large amounts of excess 
flux at high latitudes. The two disk models (solid lines) do a much 
better job fitting the vertical light distribution at all latitudes and 
show only a small systematic trend to over-subtract at high z. 



index N free, we can either fit low z or high z regions of 
the disk weU, but not both regions simultaneously; and 
(3) the two disk model is superior at fitting the vertical 
profile at all latitudes. We conclude that our galaxies 
are best modeled by the superposition of two distinct 
components with unique scale heights. 

While we already suspected the galaxies were com- 
posed of multiple stellar components based on the ob- 
served vertical color gradients and the analogous struc- 
tures present in the Milky Way, we have now shown 
that this conclusion can be derived from it!-band images 
alone. This analysis does not preclude the existence of 
additional components beyond the two disks considered 
here, however, although our data do not obviously re- 
quire them. Our fits also do not demand that the two 
components trace kinematically and chemically distinct 
stellar populations that are directly analogous to the 
thick and thin disks of the Milky Way. On the other 
hand, when combined with the color gradients observed 
in Paper II, the data are fairly suggestive of the presence 
of two genuinely distinct components. We will revisit 
this issue further in ii4.4.1l 

Ideally, we could use a statistical goodness-of-fit test 
to show that a second disk component is required when 
modeling edge-on galax ies. To establish the need for a 
thick disk in UGC 7321. lMatthew^ l)200(]|) used an F-test 
defined by 



F 



[x^(l)~X^(2)]/(p-fc) 
X'i2)/ (n-p) 



(4) 



where X^(l) characterizes the single disk model with k 
free parameters while x^(2) characterizes the more com- 
plex model with p free parameters and n total data 
points. Comparing our two disk models to single disk 
models with fixed N, the F-test favors two disks at the 
95% level or higher confidence for 32 of the 34 galaxies, 
confirming that the two disk model is a better fit than a 
single disk, as seen in Figure [7| Even if the index N is 
allowed to vary, the two disk model is still favored in 31 
of the 34 g alaxie s. Although the F-tcst works well for 
iMatthewsl (j2000() when fitting one-dimensional profiles, 
there are several caveats we must note for our sample. 



First, our models do not necessarily minimize the formal 
value because of our inverse weighting system. Sec- 
ond, the F-test relies on the formalism and thus as- 
sumes all errors are random and Gaussian. As we noted 
in t l2.2l our residuals are definitely non-Gaussian, and 
therefore the results of any F-test should be viewed as 
suggestive but far from conclusive. 

4.2. Why Not a Halo? 

In addition to modeling the galaxies as a superposition 
of thick and thin disks, we investigated models composed 
of a singl e disk and a "stellar halo" component as advo- 
cated bv iZibetti et alJ l|2004D . For our halo model, we 
used a generalized Hubble density distribution (,Wu ct all 
|2002() with the luminosity density 



-Chalolf, z) 



Co. 



halo 



{l + [r2 + (z/g)2]/r2}7/2 



(5) 



Viewed edge-on, this density distribution projects to the 
surface brightness profile 



-{z/q)T- 
(6) 

where F is the standard gamma function 

We find our data strongly prefers a second disk com- 
ponent to a halo. Over half of our halo fits converged 
on very flattened halos {q < 0.45), essentially reproduc- 
ing the properties of a thick disk, although one with a 
radial gradient in scale height. In addition, 40% of the 
fits converged to halo luminosities that are less than 1% 
of the disk luminosity, implying that the fitting proce- 
dure cannot actually use the new component to improve 
upon the single disk fits. When unconstrained, the halo 
exponential parameter 7 ran away to very large or small 
values (producing a uniform background or a compact 
point source), again implying that a power-law halo is 
not the appropriate model for the light distribution at 
high z. 

Using the F-test defined in H4.ll to compare the two 
disk fits to 9 free parameters with the disk plus halo 
fits with 10 free parameters, we found only 11 of the 
galaxies were better fit with a halo than the second disk. 
Even when suggested that the halo model was a bet- 
ter fit, the flattening parameter converged to extreme 
values (less than 0.4, or greater than 1) in 7 of the 11 
cases, thus flattening the halo into a more disk-like struc- 
ture. In those cases, the preference of a halo component 
may indicate the presence of a radial gradient in disk 
scale height. In the few cases when a preferred halo flt 
remained roughly circular, it was because the halo col- 
lapsed to flt a small central bulge or star-forming region. 
Because these regions are bright, they can greatly affect 
the formal value of x^ and the F-test will prefer the halo 
model despite no real improvement in fitting the flux at 
large scale heights. 

The poor results of our attempted halo fitting do not 
explicitly rule out the presence of a halo at lower sur- 
face brightnesses than we can detect in our images. In- 
deed, a stellar halo like the MW's would only start 
to domin ate the thick disk co mponent at iin ~ 27.5 
mag/D" ({Morrison et al.l IT997I) around z ^ IOzq (our 
fits extend to only z ^ 3 — 4zo). In M31 the stel- 
lar halo population dominates at a projected radius of 
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Fig. 8. — Tully Fisher relation derived from the single disk fits. 
Open circles show points uncorrected for internal extinction while 
solid circles have been corrected for internal extinction using the 
method of TuUv ct al. 11998). Solid lines on the left and middle 
panel show TF relations from the Hubble Key Project ISakai et al.l 
[2000 ). The s olid li ne in the right panel shows the K' TF-relation 
of iVerheiicni 120011) . The luminosities for the two disk fits show a 
comparable offset. 



^ 30 kp c and a surface brightnes s level of fxy ^ 31 
mag/D" l)Guhathakurta et al.ll200,^ . and thus comprises 
< 5% of the total stellar luminosity. As before, if a 
comparable halo component was present in our sample 
it would be much too faint to be detected in our data. 

'Zib etti et al.l l|2004D fit a composite galaxy created by 
stacking over 1000 edge-on spirals from the Sloan Digital 
Sky Survey. Using equation EJ they found a slightly flat- 
tened halo with q = 0.50 in 5, 0.60 in r and i, and 0.70 in 
z. There are several reasons we believe the excess light 
we detect at high latitudes is not equivalen t to the halo 
component discussed bv .Zibetti et alJ l)2004|) : (1) Our R- 
band is close to r and i, yet when we try to fit a halo 
component, our values o f q are much lower with a me- 
dian value of 0.4; (2) The lZibetti et all ©04) halos only 
begin to dominate the surface brightness at very large 
heights (z = IBzn.) beyond our 1-cr cropping limit; (3) 
iZibetti et a l. (2004) find that their stellar halo becomes 
prominent at a surface brightness of fXr ~ 27 mag/D", 
fainter than what we can detect in our individual images. 

4.3. Dust effects 

We have previously found that galaxies in our sam- 
ple with rotational velocities greater than 120 km s^^ 
host concentrated dust lanes ( Dalcanton et al. 2004). We 
therefore need to consider the effect that dust extinc- 
tion will have on our fitted parameters. To quantify the 
amount of extinction in our edge-on sample, we compare 
the total luminosities of our best fit models to the Hub- 
ble Key Project TuUy-Fisher relation ( Sakai ct al. 2000) 
in Figure |S| We find that all of our galaxies, even those 
without recognizable dust lanes, lie significantly below 
the TF-relation for fa ce-on spirals . If w e apply the ex- 
tinction correction of iTullv et al.l (|1998() . however, our 
extinction corrected total luminosities move nicely onto 
the facc-on TF relation. 

Their offset from the TuUy-Fisher relation implies that 
our models do not capture all the stellar flux from our 
galaxies. There are several possible ways dust could 
influence our fitted parameters to yield lower than ex- 
pected total luminosities: (1) the peak surface bright- 
nesses could be too low; (2) the scale lengths could be 
too short; (3) the scale heights could be too small; (4) 
the vertical profile could appear less peaked than it truly 
is (e.g., a sech^ instead of an exponential); or (5) some 
combination of the above. 

We can say with some certainty that the scale lengths 
do not appear to be shortened by dust extinction. If any- 



thing. Figure |2| shows that our scale lengths are larger 
than those measured in face-on systems. In a similar 
fashion, it is unlikely that our scale heights are short- 
ened greatly due to dust, as their bias is likely to have the 
same sign as the scale lengths. Moreover, our weighting 
scheme de-emphasizes the midplane region, and our scale 
height fit is therefore dominated by fiux coming from 
high galactic latitudes. Therefore, any dust distribution 
which is concentrated along the midplane, or uniformly 
distributed though the galaxy, should have little to no 
impact on our fitted value for the scale height. Only a 
truly pathological dust distribution, such as one having 
large amounts of dust at high z compared to the mid- 
plane, would skew our scale height parameter to lower 
values. 

Having eliminated biases in the scale height and scale 
length, we find that dust extinction is most likely affect- 
ing our choice of vertical profile and/or the fitted value 
of the peak surface brightness. Unfortunately, there is 
a degeneracy between these two parameters which could 
only be broken if we knew the intrinsic dust distributions 
in our galaxies. If the dust affects only the midplane re- 
gion, then the error could be confined to just the vertical 
distribution, while a more diffuse and vertically uniform 
dust distribution would lower the central surface bright- 
ness but not affect the vertical profile. 

We conclude that while our galaxies display clear signs 
of internal extinction caused by dust, the lost flux will 
cause us to either pick the wrong vertical proflle (which is 
not of particular importance since we are not concerned 
with the midplane behavior), and/or underestimate the 
overall flux normalization as parameterized by the edge- 
on central surface brightness. However, since th e em- 
pirical extinction correction of ITullv et al.l (|1998) does 
an excellent job moving our galaxies onto the TF rela- 
tion (despite the correction originally not being intended 
for use on galaxies with extreme inclination angles), we 
chose to apply this correction to the luminosity of our 
thin disk component. We do not assume any correction 
for the extended thick component, since a much smaller 
fraction of its projected area would be obscured by dust. 

4.4. Results of Thick + Thin Disk Fits 

We now discuss the results of fitting two disk compo- 
nents. We fit a total of six 2-disk models, each with differ- 
ent combinations of = 1 and N = 2 vertical profiles for 
the thick and thin components. We also considered mod- 
els convolved with a cr = 1" circular Gaussian (to model 
seeing) and models where the midplane (±zo,singie) is 
masked (to avoid dust lane contamination). The prop- 
erties of the models are described in Table 4. We use 
the inverse-weighting scheme for all the fits, as we found 
that one disk component always collapses to fit bright 
midplane structures if more conventional weighting is 
adopted. As discussed above, we fit the two-disk mod- 
els only to the i?-band images. Our i^Ts-band data does 
not go deep enough to reliably detect the thick disk com- 
ponent, and the i3-band suffers from dust extinction, is 
biased towards young stellar populations, and is a poor 
tracer of the faint red light expected from an old thick 
disk. 

The resulting parameters for the fits are listed in Ta- 
ble 5. For the central value of each parameter we list 
the median value of convergent N ~ 2 models. The 



Structure of Thick Disks 



9 



A MW youDg thin 
T MW old thin 
■ MW old thick 


+ ■ 







I 100 150 200 250 
Circular Velocity (km s"') 



I 100 150 200 250 
Circular Velocity (km s ') 



Fig. 9. — Scale he ights of thin and thick dis ks. Values 
for the Milky Way from lLarsen Hz Humphrevsl 12003) are plotted 
for comparison using zq = 2/1^ for an exponential vertical pro- 
file. Open symbols are used for galaxies with prominent dust 
lanes. Dashed lines show power-law fits to the data (zq thin = 

(610 P-)(tooT^)°-'° -o.thick = (1400 pc)(^5^^^)i-0 
). In galaxies that have strong dust lanes, the scale height of the 
thin disk is likely to be biased towards larger values. 



the galaxies or from 2-d fitting similar to the procedure 
used in this paper. We summarize these other results in 
Table 2. Figure ITTI indicates that our scale height ratios 
are slightly lower than those measured in other systems 
{zo,thick/za^thin ~ 3), implying that our derived thick 
disks may be ~ 25% thinner and/or our thin disks are 
thicker than those derived in other galaxies with other 
methods. However, ou r median ZQ^thick/ zo^thin is very 
similar to Nee ser et al.l |2002)'s measurement of the LSB 
galaxy ESO 342-017 , the most comparable galaxy in the 
literature to galaxies in our sample. These differences 
may indicate that the thick disks of early type galaxies 
may be proportionally thicker than those of late type 
galaxies. 



uncertainties are the full range of values to which the 
differ ent models in Table 5 converged, as discussed in 
ii2.2\ We also list the ratio of total luminosities for the 
model thick and thin disks. The range of luminosity ra- 
tios include models with disks having = 1 or = 2. 
The luminosity ratios are calculated only for flux which 
falls inside the 1-cr noise region (i.e. only the region that 
was included in the fit). The true luminosity ratios could 
therefore be different from our quoted values if the disks 
extend far beyond our detection limits. We have mea- 
sured the size of this correction by extrapolating the fits 
and find it can change the luminosity ratios by only 10% 
at most. The luminosity ratios in Table 5 do not include 
the extinction corrections derived in 

4.4.1. Scale heights of the thick and thin disks 

The scale heights of our thick and thin disks are plot- 
ted in physical units in Figure El The scale heights of 
both the thin and thick disks increase systematically 
with circular velocity. Fitting power laws to the re- 
lations, wc find zo^thin = (610 Pc)( ioo km s-i )°'^° ^^'^ 
zo,tMck = (1400 Pc)( ioo km s-i )^'° ^i*^ ^'^^ scattcrs of 
30% in both cases. In general, the scale heights of the 
two disks bracket the scale height derived for a single 
disk, as expected. 

For massive galaxies with large circular velocities (T4 > 
170 km/s), our derived value for the scale height of the 
thin disk is 2-3 times larger than the MW's thin disk. 
However, these galaxies have the most prominent dust 
lanes, which may substantially obscure our view of the 
thin disk. It is therefore likely that the scale heights 
of the thin disk may be significantly overestimated in 
these cases. The plot of zq.r/ zq.k for the single disk fits 
(Figure 01 lower right) is also consistent with this inter- 
pretation. Unfortunately, this limitation is unavoidable 
until sufficiently deep K-hand data is available. 

Figure [TUl shows the ratio of the thick to the thin disk 
scale height zo^tMck/ zo^thin- We find a mean ratio of 2.5 
with a scatter of 30%. In Figure ITTI we show our data 
along with other thick and thin disk scale heights de- 
rived from the literature. For the Milky Way, these scale 
heights are derived from star counts. For the other lit- 
erature values, the scale heights are derived either from 
fitting double exponential profiles to 1-d cuts through 
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TABLE 2 
Thick disks from the literature 



galaxy name 


Type 


Band 


Fitted 


Vc (kms-l) 


2q thin 






34 galaxies 


ba 


H 


z-a 


55-190 


1.6-5.5 


0.6-1.6 


0.07-7 


6 galaxies 


ba 


star counts 


i-Q 


67-131 


1.7-2.7 






ESO 342-017 


Scd 


R 


1-d 


127 


2.5 


1 


0.45 


IC 5249 


Sd 


R 


1-d 


110 


3 


0.6 




MW 


Sbc 


star counts 


2-d 


220 


3 


1.3 


~ 0.13 


NGC 6504 


Sb 


R 


1-d 


110^ 


3.9 




~0.4 


NGC 891 


Sb 


R 


1-d 


224 


3.5 




0.12 


NGC 4565 


Sb 


6660 A 


1-d 


244 


2.2 


1.4 




5 galaxies 


SO 


R and V 


2-d 


- 130^ 


2.6-5.3 


1.7-1.9 


0.33-1.0 


NGC4710 


SO 


R 


1-d 


147 


3.2 






NGC4762 


SO 


R 


1-d 


110 


4.6 






Simulation 








240 


4.7 




0.35 


Simulation 








150 


2.6 


0.63 


0.8 



Reference 



This study 
Seth et al. (2005) 
Neeser et a l. (2002) 
Abe et al. ?1999) 
^arsen & Humphrexs (2^^^ 
van Dokkum et al. (1994) 
^orrison et al. ( 1 991) 
Wu et al. (2002) 
Pohlcn ct a l. (2004) 
dfiG^^^^an dcr Kruit (1996J 
Je (iriis k van d cr Kruit (1996 ) 



Brook et al. 



120031 



^Estimated from TuUy-Fishcr relation 

^Dynamical information only available for 2 of the 5 galaxies 
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Fig. 10. — Ratios of the scale heights for the thick and thin 
disks. Error bars represent the full range of ratios to which dif- 
ferent models converged. Galaxies with prominent dust lanes are 
plotted as open circles. The dotted line shows the median value of 

^0, thick /^O, thin = 2.35. 
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Fig. 11. — Comparison of our scale height ratios to values drawn 
from the literature. We compare to early-type galaxies (Sb and 
earlier; red diamonds), late type galaxies (Sc and Sd; blue trian- 
gles), the MW (green square), and simulated galaxies. The range 
of values available for thepohlen et al. (2004) sample of early type 
galaxies is plotted as a single point with error bars. 



We also do not a priori know whether our thick and 
thin components are strict analogs of any particular com- 
ponent in the disk of the Milky Way, which is usually 
broken into at least three components; (1) the "young 
star-forming disk" (zq ~ 200 pc) which is dominated 
by molecular clouds, dust, and massive OB stars; (2) 
the "old thin disk" (zn ~ 600); and (3) the "thick disk" 
(zn kpc) dB ahcall fc Soneiralll980HReid fc Maiewskil 
|1993':'Bus er et al.ll9 99: Larsen & HumDhrevf^2003^'Oiha! 
EqOI; Che n et all 12001). which contains ^ 15% of the 
total disk light (Buscr ct glj Il999t IChen et alJ l200lt 
iLarsen fc Humphreys . 2003.) . 

One possibility is that our thin and thick disks might 
be analogous to the MW's young star forming disk and 
old thin disk, respectively, with no detectable analog of 
the MW thick disk. However, we do not believe the sec- 
ond component we have fit is an "old-thin" disk. The 
scale heights of our thin disks are larger than what has 
previously been measured for thin star forming layers. 
[Matthews (2000) find UGC 7321 has a "young disk" 
with zq ~ 185 pc. Similarly, IC 2531 has a young disk 

with zo ^ 134 pc (,Wainscoat et al 1989) and the MW's 

young disk has zn ^ 200 pc l|Bahcall fc SoneiralllQSCl: 
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Fig. 12. — Comparison of our results with the scale heights of dif- 
ferent stellar populations measured from resolved stars in 6 nearby 
galaxies I Seth et al. 2005). The component we have identified as 
the thin disk appears to be intermediate between the scale height 
of young Main Sequence stars and medium-age Asymptotic Giant 
Branch stars while our thick disk component is similar to the old 
Red Giant Branch populations. 



iReid fc Ma iewski 119931) : (using the conversion that at 
large scale heights zq ~ 2/12). All known young star- 
forming disks thus have Zq ~ 200 pc. The only galaxies 
in our sample that have thin disk scale heights approach- 
ing values this small are far less massive than any of the 
galaxies in the previous studies. 

Having ruled out an old thin disk, we now consider 
the possibility that our second thicker disk component 
is analogous to the MW's thick disk. We find strong 
support for this possibility from studie s of res olved stel- 
lar pop ulatio ns in similar syst ems (e.g.. lSethe t al. 200f!| 
lMouldll2005HTikhonov et al.! '2005'). In particular, a re- 
cent analys is of resolved ste llar populations in edge-on 
galaxies bv lSeth et al.l l)2005l) separates stars into young 
Main Sequence (MS), older Asymptotic Giant Branch 
(AGB), and still older Red Giant Branch (RGB) stars. 
ISeth et abl (2005) find that the younger stellar popula- 
tions have systematically smaller scale heights than the 
ancient RGB population. In Figure IT21 we compare our 
thin and thick disk scale heights w ith th e MS, AGB, and 
RGB scale heights of Seth ct al.| l|2005D . We find that 
our thin disk components have scale heights very similar 
to the young and intermediate age stellar populations of 
ISeth et al. (2005|), while our thick disk components have 
scale heights similar to, or perhaps slightly larger than, 
the old RGB populations. Figure^] supports that what 
we have labeled the thin disk hosts a young and inter- 
mediate age stellar population akin to the thin disk of 
the Milky Way while what we have labeled as the thick 
disk traces a different older and redder population, not 
an extension of the thin disk. When coupled with our 
observation of strong vertical color gradients (Paper II), 
and kinematic differences above and at the midplane, we 
believe there is compelling evidence that the second disk 
component required by our surface photometry does rep- 
resent a truly distinct stellar population. 

4.4.2. Ratio of scale lengths 

Physical values of the thick and thin disk scale lengths 
are plotted in Figure El We see a systematic in- 
crease in the radial scale lengths of both disk com- 
ponents with galaxy mass. The data are well fit 
by hR^thtn = (3.4 kpc)( ^,,Q s-i )^'^ ^'^'^ hn^thick = 
(^•9 kpc)( ^QQ km s-i )^ ° '^it^ scatters of 22% and 

29% respectively. 
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Fig. 13. — Scale lengths of thick and thin disks from the 2-disk 
fits. The dashed lines show power law fits. Open symbols are used 
for galaxies with prominent dust lanes. Dashed lines show power- 
law fits of /lij thin = 3.40(Vc/100 km s^"*^)^--^ kpc and ?ti{_thick = 
3.9(yc/100 km s-^)l O kpc. 
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Fig. 14. — Ratios of the scale lengths for the thick and thin 
disks. The horizontal line indicates where the thin and thick disk 
components have equal scale lengths. Error bars represent the 
full range of ratios to which different models converged, and are 
indicators of our systematic errors. Open symbols are used for 
galaxies with prominent dust lanes. 
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Fig. 15. — Comparison of our scale height ratios to values dr awn 
from the literature. Comparison points are the same as Figure [TTl 



In Figure E| we plot the ratio of the thick to thin 
disk scale lengths. We find that the thick disks have 
systematically larger scale lengths for all but 5 galaxies. 
Thick disks with long scale lengths are in excellent agree- 
ment with previous thick disk measurements, as shown 
in Figure 1151 where we include data from the literature 
(Table 2). In all but one measurement of physical (i.e., 
non-simulated) thick disk scale lengths, the thick disk is 
found to be slightly longer than the thin disk. 



We were initially concerned that this result could be a 
systematic result of our weighting and masking scheme. 
For example, if our galactic disks truncate at large radii 
as found in other edge-on sy stems (Krcgelet al. 200^ 
IKreeel fc van d er Kruit'20041 , then our model fits would 
converge to have the fainter thick disk dominate at large 
R. However, we included a model (Table 4) where the 
midplane is masked (which also effectively removes re- 
gions of the galaxy where disk truncation would be de- 
tectable) and still found that the thick disks have longer 
scale lengths. 

We note that there are some limitations in interpreting 
our scale lengths, particularly for the thin disk. First, 
the derived radial scale lengths do not necessarily re- 
flect the stellar radial scale length. The thin disk in 
particular shows a strong radial color gradient, imply- 
ing a mass-to-light ratio that decreases with increasing 
radius (see color maps in Paper II). This trend suggests 
that the radial scale length of the stellar mass should 
be even smaller for the thin disk, further increasing the 
ratio h}i,thick/h Ruthin- We may also have overestimated 
the scale length of the thin disk if it is affected by dust 
in a manner similar to the what is observed in our single 
disk fits (Figure 13 and . Both of these effects sug- 
gest that hii^thick/hR^thin may be even larger than indi- 
cated by Figure^] On the other hand, HI is typically 
more extended than the optical disk ijSwaters et al.l20"02 
iBeeum et al.ll2005j) . such that the radial scale length of 
the baryons in the thin disk may be longer than indicated 

by /i_R,thin- 

4.4.3. Axial ratios of the thick and thin disks 

The axial ratios {hn/ zq) for our thick and thin disks are 
plotted in Figure ^1 along with values for the MW thick 
and thin components for comparison. We find our thick 
disks have a mean Hr/zq = 3.4 with an RMS scatter of 
1.7 while the thin disk has a mean value of hfj/zo =4.7 
and RMS scatter of 1.8. Our thin disks therefore tend 
to be comparable to the MW thin disk, but are slightly 
rounder at low masses, in agreement with other studies 
(see 0. The axial ratios of the thick disks show a large 
spread in axial ratios, and are in general comparable to, 
or slightly thicker than the MW thick disk. However, 
the radial scale length of the MW is not well constrained 
since it is determined primarily from star counts near 
the solar circle. We also note that the thick disk com- 
ponent is drastically rounder than the MWs old thin 
component, further ruling out the old thin disk as an 
explanation for our second disk component. 

4.4.4. Peak surface brightnesses 

The edge-on peak surface brightnesses for our two disk 
components are plotted in Figure lTTI There is a trend for 
more massive thin disks to have brighter peaks, similar 
to the trend seen in the single disk fits (Figure El ■ The 
thick disk components show a large amount of scatter in 
their peak values. 

Performing a naive transformation to convert to the 
face-on orientation, the central surface brightness be- 
comes fio — /i(0, 0) — 2.51og(zo/^fl)- We find that the 
average central surface brightness of the thick disk is 0.6 
mag/n"fainter than the thin disk, implying only ~ 35% 
of the stellar flux in the i?-band would come from the 
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Fig. 16. — Axial ratio (hn/zo) vs circular velocity for the thin 
(left) and thick (right) disks. Galaxies with prominent dust lanes 
are plotted as open circles. The axial ratios of Milky Way disk 
components are plotted for comparison (Table 2). The axial ra- 
tios of our thin and thick disks agree well with the comparable 
components for the MW. 
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Fig. 17. — Edge-on peak surface brightnesses for the 2-disk fits 
in the _R-band. Galaxies with prominent dust lanes are plotted as 
open circles. 



thick disk if the galaxies were viewed face-on. In the more 
massive galaxies, the face-on central surface brightness of 
the thick disk can be up to 2 mag/D" fainter than the 
thin disk. These values do not include corrections made 
for extinction. Presumably, the thin disk would suffer 
less extinction when viewed face on, and would further 
dominate the observed stellar flux. It is therefore not 
surprising that the thick disk is largely undetected in 
face-on galaxies. 

The total integrated colors of the galaxies will be bi- 
ased towards the thin disk population as well. After 
making the extinction corrections in H4.3I we find the 
total integrated colors of our low mass-galaxies (14 < 
100 km s~ ) are in the range 0.5 < B — R < 1, much 
bluer t han the thick disk (see F igure EOI below). Us- 
ing the lBruzual fc CharlotI l)2003}) stellar synthesis code, 
these colors correspond to a stellar population burst with 
and age of ~ 1 Gyr, or a galaxy with a uniform star for- 
mation history. Thus, in spite of the substantial thick 
disk population, the mean colors of the galaxy reflect 
only the youngest disk population. 

4.4.5. Ratio of luminosities 

We now compare the total luminosities of the thick and 
thin disks (Figure ITS)) . In our raw fits, the luminosity of 
the thin disk is almost certainly underestimated due to 
the effects of dust, as shown in H4.3I To correct for dust, 
we assume that all flux lost from extinction (Figure IHJ 
should be assigned to the thin disk. This correction will 
give us the most conservative estimate for the contri- 
bution of the thick disk to the total stellar luminosity. 
Figure 1181 shows a strong trend with mass (Spearman 
p — —0.70, 4. Oct). Thick disks of high mass galaxies 
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Fig. 18. — Ratio of the total i?-band luminosity of the thick 
disk compared to the thin disk for the sample galaxies. The thin 
disk luminosities have been corrected for internal extinction. The 
dotted line indicates where the thick and thin disks contribute 
equally to the total luminosity. Error bars show the full range 
of values from different models, and are indicative of our sys- 
tematic errors. Galaxies with prominent dust lanes are plotted 
as open circles. The solid line is a fitted power-law of the form 
Wk/ithitr = 0.25(y,/100 km s-l)-2-i. 
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Fig. 19. — Comparison of our luminosity ratios to others in the 
literat ure. Comparison points are the same as Figure ITTl Unlike 
Figure 1181 we have made no correction of internal extinction to 
allow for easier comparison with previous studies. 



{Vc > 120 km s-i) contribute ~ 10% of the total lu- 
minosity of the galaxy, while in lower mass systems the 
thick disk contributes up to 40% of the total luminos- 
ity. This trend can be well represented by the relation 
-^thick/ithin = 0.25(14/100 km s~^)~^-^, shown as asoHd 
line in Figure [THl 

We compare our measurements to previous thick disk 
measurements in Figure lHH Unfortunately, there are few 
measurements of total disk luminosities in the literature. 
When possible, we have taken other authors' disk pa- 
rameters and calculated the resulting total luminosities 
(see Table 2). For the Milky Way, the local stellar den- 
sity of thick disk stars has consistently been measured 
between 4 and 10 % of the lo cal thin disk density (e.g., 
iBuser et al.. 1999.: Chen et al.„ 2001) . which corresponds 
to a total luminosity ratio of ~13% for reasonable esti- 
mates of scale heights, lengths, and mass-to-light ratios 
for the two disks. Because the values of -/jthick/^thin from 
the literature do not include internal extinction correc- 
tions, we compare them to our uncorrected luminosity 
ratios. 
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Figure 1191 shows that our luminosities compare well 
with other thick-thin disk systems in the literature. The 
higher mass galaxies in our sample {Vc ~ 140 — 200 km 
s~^) tend to be thin disk dominated with ithick/^thin 
0.1 — 0.2 (corrected for extinction), like the Milky Way 
and NGC 891. Intermediate mass galaxies (70 < Vc < 
100 km s~^) have slightly more luminous disks , simi lar to 
measurements of ESQ 342-017 (iNeeser et al.ll2002|) and 
SO's fPohlcn et al. 2004). Unfortunately, we cannot find 
any comparable measurements of thick disks in the low 
mass systems (50 < Vc < 70 km s^^) that are thick disk 
dominated in our sample. 

We believe our measurement of the total luminosities 
are more robust than the measures of the peak surface 
brightness . Central surface brightnesses depend strongly 
on the vertical profile and can vary greatly from au- 
thor to author. On the other hand, our fits of the to- 
tal luminosity are good matches the the data (Immodoi 
- mobservedl ~ 0.2 mags), and fall on the Tully-Fisher 
relation (Figure |HJl. 

4.4.6. Mass ratios 

Figure [T51 indicates that thick disk stars provide a sig- 
nificant fraction of a galaxy's total luminosity. However, 
as seen in Paper II, the thick disk tends to have a redder 
color than the thin disk, especially in low mass galax- 
ies, and thus will have larger stellar mass-to-light ra- 
tios than the thin disk. Therefore, the thick disk may 
well dominate the stellar mass in many of our galaxies. 
We estimated the stellar disks' masses using the lumi- 
nosities of the two disk components, along with color 
information from our single disk fits. Specifically we 
used the spectrophoto metric galaxy evolution analysis of 
IBell fc de Jond l|200H) to convert our B — R color maps 
into stellar mass-to-light ratios for each disk, and then 
convert our luminosity ratios into mass ratios for the 
thick and thin components. 

The initial analysis of vertical color gradients in our 
sample (Paper II) suggested that the colors of thin disks 
vary systematically with galaxy mass, but that the colors 
of thick disks are fairly uniform. We therefore assumed 
that the thick disks have uniform colors and mass-to- 
light ratios in each galaxy. To convert disk colors to 
masses, we first analyzed our it!-band 2-disk fits to find 
regions where the thick disk contributes more than 75% 
of the total flux inside the l-cr noise contour. Out of 34 
galaxies, 27 have a clearly thick disk dominated region. 
We created a model B — R color map of each galaxy 
using our 2-disk J?-band model and single-disk i3-band 
fits. Using the model images allowed us to avoid dust 
lanes, HII regions, and foreground objects that would 
complicate an analysis on the real images. We then used 
this color map to find the average B—R value in the thick 
disk dominated regions and took that as the approximate 
color for all the thick disk stars. We assumed the thick 
disk has a constant color, thereby guaranteeing that its 
structural parameters will be the same in both the B and 
R bands. With the B — R thick disk color from the model 
color map and the thick disk structural parameters from 
the 2-disk fit, we then made a model B — R color map for 
the thin disk by subtracting off the thick disk component 
from both the _B-band and _R-band models. 

We applied internal extinction and reddening correc- 
tions to our models using the results of H4.3l and assumed 
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Fig. 20. — Extracted colors for the thick and thin disks. The thin 
disk colors shows the full range of B — R values for the midplane 
between hn < R < 3hfi. The thin disk has been corrected for 
internal extinction, but we assume no correction for the thick disk. 
Open circles show galaxies with dust lanes. 




Fig. 21. — Comparison of thin and thick disk colors. The di- 
agonal line indicates where the two components have equal color. 
Thin disk colors have been corrected for internal extinction. Open 
circles show galaxies with dust lanes. 



that dust had a uniform effect on the thin disk colors but 
a negligible effect on the thick disk. Using this approx- 
imation, we found E{B — R) ^ 0.1 for low mass galax- 
ies and ~ 0.4 for higher mass galaxies. Although it is 
only a rough approximation, our reddening correction is 
in good agreement with t he radiative transfer model of 
iMatthews fc WoodI l|200H) who find that most disk light 
in their modeled edge-on galaxies suffer reddening of or- 
der E{B — R) ~ 0.1, and that the reddening saturates at 
E{B-R) = 0.31. 

The resulting colors for thick and thin disks are plot- 
ted in Figure [201 The thick disks tend to be red with 
1.0 < B — R < 1.7, while the thin disks are blue in 
low mass galaxies and become nearly as red as the thick 
disks in the higher mass galaxies. This trend is also seen 
in Figure 1211 where we directly compare the colors of 
each component. 

Using the colors shown in Figure EOl a nd the color de- 
pende nt stellar mass-to-light ratios from IBell fc de .Tonal 
()2001|) . we converted the thick and thin disk luminosi- 
ties to stellar masses using M — {M/ L)]iLfi, where 
Lfl is the extinction corrected i?-band luminosity from 
^4.4.51 For the thin dis k, we calculated {M/L)ji using the 
iBeff fc de Jo^ l|2001fl model which assumes a Salpeter 
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Fig. 22. — Stellar mass ratios of the thick and thin disks. The 
thin disks luminosities and corresponding masses have been in- 
creased to account for dust extinction. Galaxies with prominent 
dust lanes are plotted as open circles. The solid line is a power-law 

fit (Mthick/A/thin)stars = 0.53(^/100 km S-l)-2-3. 
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Fig. 23. — Baryonic mass ratios of the thick and thin disks. 
The thin disks luminosities, and corresponding masses, have been 
increased to account for extinction effects and include the esti- 
mated thin disk mass stored in gas. Galaxies with prominent dust 
lanes are plotted as open circles. The solid line is a power-law fit 

Mthick,stars/MtMn,baryons = 0.19(^/100 km S"!)"!'*. 



IMF and metallicity of Z — 0.02, for the thick disk we 
use the same model with Z = 0.08. Overall, our results 
were insensitive to the s tellar evolution and m etallicity 
differences covered in the lBell fc de Jon3 II2001D models. 

The resulting mass ratios of the thick and thin disks 
are plotted in Figure!^ As expected, Figurel^confirms 
the features from our luminosity analysis. First, there is 
a strong trend for lower mass galaxies to have a larger 
fraction of their stellar mass in a thick component. The 
trend has a Spearman correlation of p = —0.86 (4.1(t) 
and can be well fit with the relation Mthick/A^thin = 
0.53(14/100 km s"^)"^-^. Second, in low mass galaxies, 
~ 1/3 to greater than 1/2 of the stellar mass is in the 
thick disk. Thus, the stellar mass of very low mass galax- 
ies are dominated by thick disk stars. 

Part of the trend in Figure |221 may be due to low star 
formation efficiency in lower mass disks. These systems 
have high gas mass fractions, and thus may not yet have 
built up a significant stellar mass in the thin disk. To in- 
vestigate this possibility, we calculated the baryonic mass 
fraction of the thick and thin disks, assuming that all gas 
in the galaxies is associated with the thin disk and that 
the thick disk is entirely stellar. We calculate the gas 
mass as Mhi/Mq = 236(P J SdV where d is the distance 
to the source in Mpc, and S is the flux density i n mJ y 
over the profile width dV in km s^^ ( Zw aan et aliil997fl . 
To account for He and metals, we make the standard cor- 
rection Mgas = 1.4Mhi- We do not include a correction 
for molecular gas. 

Figure 1231 shows the resulting baryonic mass ratio of 
thick and thin disks with the mass of HI gas added to 
the thin disk component. When the gas is included in the 
thin disk component, we find that none of the galaxies 
remain thick disk dominated although the baryon mass 
fraction in the thick disk does remain substantial for low 
mass galaxies. Eleven of our galaxies had no HI data 
and their gas fraction was estimated by fitting a sim- 
ple power law to the gas fraction of our galaxies with 
HI measurements. Figure |^ shows the final calculated 
baryon fractions for all of the stellar and gaseous com- 
ponents, and clearly indicates the increasing importance 
of the thick disk in lower mass galaxies. 
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Fig. 24. — Final baryon fractions of the thick and thin compo- 
nents as a function of maximum velocity. The thin disk has been 
corrected for internal dust extinction. The baryonic mass fraction 
of the thick disk clearly increases towards lower galaxy mass. 



5. THE FORMATION OF THE THICK DISK 

Given evidence from the Milky Way and from nearby 
resolved galaxies, we will assume in the following discus- 
sion that thick disks have a formation mechanism dis- 
tinct from that which forms the thin disk. We will also 
assume that the properties of "thick" and "thin" com- 
ponents from our 2-D fits will roughly approximate the 
properties of the chemically and kinematically distinct 
thick and thin disk analogs of the Milky Way. As we 
have argued above in H4.4.1I the structural properties 
of the fits are consistent with those of the corresponding 
components of the Milky Way and with the results of de- 
tailed stellar population studies in nearby galaxies (e.g. 
Figure I12|) . We therefore will simply assume a perfect 
correspondence between our fits and distinct thick and 
thin components for the rest of this discussion. While 
this assumption is not ideal, it is unavoidable, given that 
a full kinematic and chemical analysis of the stellar com- 
ponents is essentially impossible far outside the Local 
Group. 

5.1. The Merger Origin of Thick Disks 
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As discussed above in 21 there are three general classes 
of thick disk formation scenarios - one where the thick 
disks stars form in situ, one where the thick disk stars 
form in a thin disk but are then impulsively "heated" 
to large scale heights, and one where thick disk stars 
form first in other galaxies but are then directly ac- 
creted. In the last decade there has been a growing 
body of evidence in favor of the latter two scenarios. 
This evidence includes the detection of strong kinematic 
diff erences between thic k and thin disks in other galax- 
ies (lYoac him fc Dalcanton 2005) and in the Milky Way 
(iCrilmore et al.ll20n2D . as well as evidence from chemical 
abundance studies for extended star formation histories 
of Milky Way thick disk stars (e.g.,^cnsbv ct al. 2005). 
These latter two scenarios are also naturally accommo- 
dated in current theories of hierarchical structure forma- 
tion, where mergers and accretion are common. 

Of the two merger-driven scenarios, we believe that the 
data favor an accretion origin for thick disk stars. The 
strongest evidence comes from our previous measure- 
ments of thick and thin disk kinematics in two late-type 
disks. We found that thick disk stars are rotating with 
only a small fraction of the rotational velocity of thin disk 
stars, and are indeed counter-rotating in one of the two 
cases studied. In contrast to the observed behavior, sim- 
ulations show that disks heated by satellites would have 
nearly the same angular momentum as the initial thin 
disk. Therefore, barring the unlikely possibility that the 
thin disk reforms from subsequent accretion of gas with 
angular momentum opposite to the original disk, cre- 
ating counter-rotating or slowly-rotating thick disks via 
vertical heating alone is problematic. Additional sup- 
porting evidence comes from recent H ST studies of re- 
solved stars in nearby edge-on gal a xies JSeth et al.ll2005l: 
IMouldll200l . While ISetlTetan II2005D finds evidence 
for some steady vertical heating, the oldest population 
of RGB stars shows no evidence of the vertical metal- 
licity gradient expected if they were dominated by stars 
that had participated in the steady heating. Instead, 
the thick population of RGB stars must have been es- 
tablished early, when the merger rate was still extremely 
high. At these early times, the hypothetical picture of 
a well-defined stable thin disk impulsively heated by a 
single merging event seems inappropriate. In summary, 
while the evidence does not yet conclusively rule out a 
single, thin disk heating merger event as the origin of 
thick disk stars, we consider it to be sufficiently unlikely 
that we will focus on interpreting the results in this paper 
in terms of the accretion scenario. 

The idea that stars well above the Galaxy's midplane 
may have be en direc t lv acc reted from satellites was first 
discussed by 'Statler (1988)'*, and then codified as a dis- 
tinct formation mechanism for t he thick disk s hortly 
thereafter in the revie w article bv lGilmo re et alJ ljl989f) 
and as "Model 7" in IMaiewskil l|l993|) 's review. Re- 
cent detailed studies of stellar structures in the Milky 
Way and M31 find evidence that stars are re gularly ac- 
creted by massive galaxies. For example, iMartin et al.l 
()2004|) find asymmetries in the distribution of M-giant 

Statler's work refers to explaining the kinematics of "halo" 
stars, but given the current understanding of Galactic structure, 
the specific kinematic data he was trying to explain were measure- 
ments of what we would now call thick disk stars. 



stars (e.g., the G alactic 'Ring' iNewberg et alJ 120021: 
lYannv et alJl200'3l) that are w ell explained bv a s ingle 
dwarf galaxy accretion event. IMartin et al.l l)2004|) also 
note that their modeled accretion event results in ac- 
creted stars having orbits similar to thick disk stars, and 
that the thick disk may be continually growing through 
in-plane accretion of dwarf galaxies. M31 also appears 
to also be actively disr upting dwa rf galaxies with an ob- 
servable stell ar stream ijlbata et a l. 2004) and a large ex- 
tended disk l|Ibata et alJl2005|) . 

Direct accretion as the dominant origin of thick disk 
stars has recently been revived by several numerical stud- 
ies of disk galaxy formation. Analyzing an N-body sim- 
ulation of a moderat ely low mas s spira l formed in a 
cosmological context, lAbadi et al.l (|2003) found a well- 
populated thick disk, more than half of which was made 
up of stars originally formed in accreted satellite galaxies. 
Although by no means definitive due to the simulated 
galaxy's unrealistically large bulge. lAbadi et a l. (2003)'s 
work lead to a revival of the notion that thick disk stars 
may have formed outside their host galaxies. Subsequent 
simulations of a collapsing sphere seeded with perturba- 
tions by Brook et al. (2004, 2005) generated thick disks 
ass ociated with an earl y episode of chaotic m erging. Un- 
like lAbadTeiai] l)2003|) . lBrook et alJ l)2004|) argued that 
the thick disk stars formed in situ from large velocity 
dispersion gas deposited by the satellites as they merged 
together to the final disk structure. Both simulations 
have significant limitations, making it impossible to de- 
cide in favor of either scenario at this time, but both 
stress the importance of merging and accretion in set- 
ting properties of the thick disk. 

In the context of hierarchical galaxy assembly, the 
above simulations point to a straightforward picture 
of disk formation that necessarily leads to the forma- 
tion of thick disks. At high redshift, galaxies exist 
largely as a collection of sub-galactic fragments. These 
fragments consist of gravitationally bound dark mat- 
ter "mini-halos" , many of which presumably host some 
amount of baryonic material. Because these systems are 
high in the merging hierarchy, they would be expected 
to be relatively dense, and thus some of the gas hosted 
by these sub-units is likely to form stars. Early on, the 
merging rate will be very high, and as these fragments 
come together, their orbits will tend to cir cularize , align , 
and decay due to dynam ical friction, as in jStatleil IpM) 
and early si mulations bv lQuinn & Goodm an (19861 and 
I Walker et al... 1,1996.) . The merged subunits will form a ro- 
tating flattened structure provided that the net angular 
momenta of the satellites is sufficientl y high. When th e 
merging rate declines sharply (z ^ 3) l|Zhao et al.ll2003l ). 
the disk will be left in place as a long-lived structure 
relatively unperturbed by significant accretion events. 
Any dense gas associated with the pre-galactic fragments 
must then either form stars in a b urst during the final 
merger of the fragments, as in the iBrook et a l. ( 2004 
simulations, or cool into a thin disk which later converts 
into stars. 

Within this scenario, any stars that formed in the sub- 
units and were not tidally stripped at large radii must 
necessarily wind up in a thickened disk structure, with a 
vertical velocity dispersion equal to or greater than the 
velocity dispersion of the typical pre-galactic fragment. 
Because they are effectively collisionless, the accreted 
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stars cannot lose energy and cool into a thinner disk, 
and must retain a large fraction of the initial velocity dis- 
persion and angular momentum of the satellite in which 
they formed. With this in mind, it seems impossible to 
imagine not forming a thick disk (unless star formation 
was completely supp ressed at early time s, for exampl e 
by reionization, e.g., iBullock et aI]l200(]HGnedinll200(]|) . 
The only other possible destination for the accreted stars 
would be the bulge or stellar halo. However, the sample 
considered here is essentially bulgeless. We also find no 
evidence for a luminous stellar halo down to our limiting 
surface brightness. Taking a conservative estimate for 
the surface brightness of the brightest stellar halo that 
could be present, but undetected in our data, we find 
that any stellar halo must be less than 15% of the lumi- 
nosity of the thick disk. This estimate suggests that the 
majority of directly accreted stars settle into the thick 
disk. 

In addition to the theoretical arguments for form- 
ing thick disks via direct accretion of stars, there is 
a slowly growing body of observational evidence for 
this process seen in situ at high redshift. First is 
the an alysis of high redshift " clump -cluster" galax- 
ies by lElmegreen fc ElmegreenI l)2005ft . Morphologi- 
cally, these galaxies appear to consist of many dis- 
tinct, hig h surface brigh t ness c lumps merging together. 
lElmcgrccn fc ElmegreenI (|2005D argue persuasively that 
these systems will wind up in a thickened disk with high 
velocity dispersion, and are thus likely precursors to thick 
disks. The colors of the clumps suggest that they al- 
ready contain some stars , and a re not pure gas systems. 
lElmegreen fc ElmegreenI l|2005j) also find field counter- 
parts of the clumps, suggesting that some of the stars 
may have formed before being accreted into the galaxy. 
The se cond piece of evidence is the kinematic study of 
lErb et al. ( 2004), who find that lumpy disk- like struc- 
tures at z ~ 2 show little net rotation. If their sample con- 
sists pri marily of the edge-on counterpa rts of the galax- 
ies in the lElmegreen fc ElmegreenI ()2005|) study, then the 
lack of strong rotation would be consistent with what is 
expected for material that forms a thick disk. Although 
this study traces Ha kinematics only, and thus leaves the 
kinematic state of any associated stars unconstrained, it 
would be peculiar if any stars associated with the accret- 
ing gas did not show similarly perturbed kinematics. 

5.2. Constraints from the Structures of Thick and Thin 

Disks 

In the above accretion scenario, the properties of the 
thick disk are fixed primarily by the stellar content and 
orbital properties of the pre-galactic fragments which 
merge to form the final stable disk. The properties of 
the thin and thick disks are then set by the kinemat- 
ics and gas mass fractions of the pre-galactic fragments 
when they merge. 

Within this scenario, we now discuss the implications 
of three significant properties of thick disks uncovered by 
our data: first, that thick disks are a ubiquitous and nec- 
essary component in modeling late- type edge-on galaxies; 
second, that the stellar mass of the thick disk is increas- 
ingly dominant in lower mass galaxies; and third, that 
thick disks have systematically larger radial scale lengths 
than thin disks. 



5.2.1. The Ubiquity of Thick Disks 

The 2-D fits of our sample confirm the initial sugges- 
tion of Dal canton fc Bernstei n ( 2002) that thick disks are 
a ubiquitous component of disk galaxies. Essentially all 
(32 of 34) of the galaxies which were suitable for fitting 
were statistically significantly better fit by a second disk 
component (e.g. Figure [TJ. This result adds to the pub- 
lished detections of thick disks in earlier type systems 
(summarized in Table 2). Thick disks are now routinely 
discovered in every galaxy that has been searched for 
them^ . 

The evidence therefore supports the idea that thick 
disks are a generic property of all galaxies with disks, 
from SO's to Sm's, from high masses (T4 250 km s^^) 
to low {Vc 50 km s~ ). Thick disks must therefore 
be a natural by-product of disk galaxy formation, inde- 
pendent of the formation of a bulge. The ubiquity of 
thick disks can be easily explained if most thick disk 
stars are directly accreted from pre-galactic fragments. 
As we argue above, if any star formation has taken place 
in the fragments, some fraction of those stars must wind 
up in a thick disk. The only way to avoid depositing 
the stars in a thick disk would be if the fragments were 
completely tidally disrupted at large distances from the 
central galaxies. However, at large distances the matter 
density should be much lower than in the dense cores 
of the low mass galactic fragments, making it unlikely 
that every merging satellite would experience complete 
disruption. 

The existence of widespread thick disks also suggests 
that there has been ample star formation in the very low 
mass halos which merge together to form larger galaxies. 
Most sub-units must have established stellar populations 
before merging. If instead the sub-units were entirely 
gaseous, disk galaxies would have only a thin disk com- 
ponent. Thus, there cannot have been total suppression 
of star formation by reionization up until the epoch of 
thick disk formation. 

Finally, the pervasiveness of thick disks also presents 
an additional problem for merger heating scenarios. It 
seems unlikely that every galaxy in our sample would 
have had both a merger that created a thick disk and 
accretion that reformed a thin disk. If merger heating 
was the primary driver of thick disk formation we would 
expect to find some galaxies that were able to avoid a 
destructive merger, or that failed to subsequently reform 
a thin disk. Instead, all of our galaxies require both thin 
and thick disk components. 

5.2.2. The Increasing Importance of Thick Disks in 
Lower Mass Galaxies 

In the merging picture we have adopted, sub-galactic 
fragments contribute both stars and gas to the final 
galaxy. The stars wind up in the thick disk, and the 
gas settles into the thin disk, where it gradually converts 

^ The one exception is NGC 4244, which iFWeTail 11993^ ana- 
lyzed using fits to 1-D cuts of the vertical ij-band light distribution. 
Based on th e lack of a clear break in the vertical surface bright- 
ness profile, IFrv et all 1199EI) claimed there was no thick disk in 
this galaxy. However, subs equent analyses o f the resolved stellar 
population in NGC 4244 bv lSeth et ahl <2(X)R|) and lTikhonov et alJ 
12005") revealed the presence of a clear extra-planar population 
dominated by old red giant branch stars, whose global distribu- 
tion was characteristic of a thick disk. 
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into stars. From Figure |^ we see that low mass disk 
galaxies have roughly 25% of their baryonic mass locked 
up into thick disk stars, while massive galaxies have only 
10%. Figure therefore implies a systematic variation 
in the gas richness of sub-galactic fragments at the time 
disks coalesce. In massive late-type galaxies, 90% of the 
baryonic mass must have remained gaseous during disk 
assembly, while in low mass galaxies only 75% had not 
yet converted to stars. 

Note that while we are stressing the accretion of stellar 
material to form the thick disk, our results prove that the 
vast majority (75-90%) of baryonic accretion must have 
been gaseous. If some fraction of thick disk stars did 
form in situ as suggested by Brook, then the fraction of 
gaseous accretion must have been even higher. 

There are several ways that lower gas mass fractions 
in the precursors of low mass galaxies may be achieved. 
One possibility is that the transformation of gas into 
stars proceeded further by the time the low mass disk 
galaxies coalesced. This more complete transformation 
in low mass disks could be due either to a later epoch 
of assembly, or to higher gas densities and thus higher 
star formation rates in the precursor clumps. However, 
in a closed box model, the resulting thick disk stars 
would have higher metallicities. In contrast, the esti- 
mates of the met allicities of extra-planar, RGB stars in 
iSeth et all l)2005D suggest that the metallicities of the 
thick disk stars are systematically lower in lower mass 
galaxies, compared to the Milky Way. We therefore rule 
out the possibility that star formation was more "com- 
plete" in the precursors of lower mass galaxies. 

Supernova feedback is an alternative pathway to the 
preponderance of less gas rich sub-units in low mass 
galaxies. Much of the disk material was initially in sev- 
eral subunits that were necessarily of lower mass than the 
final galaxy. Thus, the merging fragments must have had 
lower escape velocities, allowing supernova-driven winds 
to more effectively drive gas and metals from the galaxy 
at this early stage. The increased efficiency of SN winds 
in the sub-units would simultaneously decrease the gas 
mass fractions and maintain low metallicities in thick 
disk stars in low mass galaxies^. 

We can estimate the amount of gas loss needed to pro- 
duce the observed trends as follows. First, we assume 
that the observed baryon fraction in the stellar thin disk 
and gas component of massive galaxies 90%, Fig- 
ure 124(1 is indicative of the gas to stellar mass fraction 
in subgalactic fragments that are too massive to experi- 
ence significant SN blowout. We then assume that the 
precursors of lower mass galaxies lose enough gas to bring 
their gas to stellar mass fraction down to ~ 75% at the 
time of disk assembly. These simple assumptions imply 
that the sub-units of low mass galaxies must have lost 
60% of their initial baryonic mass. 

Because we have ignored possible tidal stripping of 
stars during galaxy assembly, the actual amount of gas 

® Note, however, that the overall gas mass fraction of low mass 
galaxies can remain high to the present day. The disks of lower 
mass galaxies have systematically lower baryon surface densities 
(e.g., ISwaters erallfeooa IHunter Hz ElmegreenI l200ll . and thus 
are inefficient at converting gas into stars. Their low star formation 
rate thus allows them to have higher gas mass fractions today, even 
though they were comparatively gas poor at the time their disks 
were assembled. 



lost from the precursors of low mass disks may be dif- 
ferent from what we have estimated above. However, 
assuming that tidally stripped stars wind up in a stellar 
halo, we expect the total stellar mass lost to stripping to 
be small. The Milky Way's thick disk contains a factor of 
~ 10 times more stars than its stellar halo, and thus, any 
correction due to tidal stripping is likely to be negligible. 

While our mass-dependent blowout scenario explains 
our data well, it is not clear that pre-galactic frag- 
ments actually suffer ~ 60% baryon losses due to SN 
winds. There are a wide range of results on how ef- 
fective SN winds should be at driving baryon outflow. 
At the one extreme, several groups argue that large 
SN drive n outflows ex i st in all galaxies with Vr < 100 
km s-i lIDekel fc Silkl 1198^ IPekel & Wool l20(B . At 
the the other extreme, simulations flnd that galax- 
ies with M > lO^Mf?^ e xperience almost no outflow 
l|Mac Low fc Ferrarall999j) . Similarly, o bservational con- 
strain ts on the extent of outflow vary. iMaver fc Moord 
use the baryonic TuUy-Fisher relation to claim 
that dwarf galaxies do not su ffer large removal of baryons 
while IsTrickland et alJ l|2004|) observe x-ray halos around 
massive star forming galaxies (M - 10^° - lO^M©) 
which suggest they must have ejected at least some ma- 
terial. Because we are considering the role of blowout in 
low mass progenitors of our galaxy sample, we claim that 
the current knowledge of gas blowout is moderately con- 
sistent with our scenario and we await a more definitive 
cosmological simulation which incorporates star forma- 
tion and feedback for detailed comparison to our model 
(Stinson et al., in prep.). 

There are several limitations with the simplified anal- 
ysis we have presented above. First, we have ignored 
the difficult question of how much material is accreted 
in continuous cold flows rather than bound in halos 
(Birnboim & Dckcl 2003; Keres ct al. 2004). Cold accre- 
tion of gas will tend to increase the baryonic fraction of 
the thin disk. Neglecting steady gas accretion therefore 
leads us to overestimate the gas richness of the merging 
pre-galactic fragments. Second, we have not explicitly 
considered how bulges are formed in the scenario dis- 
cussed in ii5.ll but we presume it involves repeated merg- 
ers of gas rich sub-units with little net angular momen- 
tum, or a higher frequency of major mergers in higher 
mass galaxies. Within the sample we have studies here, 
this omission is acceptable. However, more theoretical 
and observational work must be done to understand the 
thick disk population in earlier type galaxies^ 

Fina lly, we flnd it difficult to reconcile the lBrook et alJ 
(|2004') formation scenario with the increasing fraction 
of thick disk stars in lower mass galaxies. Brook ct aO 
|2004) suggest that thick disk stars form in situ from 
high velocity dispersion gas during the coalescence of 
sub-galactic fragments into a final disk. However, we 
see no obvious mechanism that could lead this scenario 
to produce a larger fraction of thick disk stars in low 
mass galaxies. One would need to invoke a mechanism 
to increase the efficiency of star formation at lower galaxy 
masses during mergers, while keeping star formation inef- 
ficient at later times. An alternative solution would be if 
steady cold flow gas accretion is more important in mas- 
sive galaxies. However, massive galaxies are more likely 
to have established a hot shock-heated halo th at would 
block cold flow (e.s.. iDekel fc BirnboinJ 120041 ^ . Thus, 
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the hkely behavior of cold flow accretion has the oppo- 
site sign as what is needed to explain the high baryonic 
fraction of thick disks in low mass galaxies. Further sim- 
ulations will help constrain this and other possible solu- 
tions. 

5.2.3. The Scale Lengths of Thick & Thin Disks 

Our data contribute to a growing number of ob- 
servations finding that thick disks have larger scale 
lengths than their embedded thin disks (Oiha 20011 
I Larsen fc Humphrevs 2003: ,Wu et al. 2002: Pohlen et all 
12004 see our Figure lOl and Table 2). The large scale 
lengths of thick disks argue against their being formed 
via vertical heating of a thin disk. N-body simulations 
find that while minor mergers can vertically heat a disk, 
they do not increase its scale length ( Ouinn et al. 199>^. 
Such mergers also tei id to leave the galax y looking like an 
earlier Hubble type l) Walker et alJll99'^ while all of our 
galaxies have no promin ent bulge cornponeii ts. As an 
example, simulations by lAguerri et alJ 1)200 If ) find that 
minor mergers can extend the scale length of the thin 
disk somewhat, by 10-60%. However, the same simula- 
tions also produce a large bulge, which is incompatible 
with our sample. 

In contrast, in the accretion scenario one would expect 
the scale length of the thin disk to be somewhat smaller 
than that of the thick disk. If the thin disk forms later 
from gas which has contracted further into halo than the 
thick disk stars, it should have a smaller scale length. If 
angular momentum is largely conserved, then the thin 
disk should also be rotating somewhat faster than the 
thick disk because of its extra contraction. 

The satellite accretion model therefore suggests that 
there may be correlations between the radial scale 
lengths and the kinematics of the thick and thin disks. 
Results in Section I4.4.2 | and photometric decornposi- 
tions by o thers (Oiha 200lt ILars en fc Humphrevsll2003t 
iWu et alJ I200a: .Pohlen et all l2004() suggest that scale 
lengths of thick disks are roughly 30% longer than those 
of their embedded thin disks, on average. Simple an- 
gular momentum conservation would then suggest that 
the thick disk should rotate with approximately 2/3 the 
speed of the thin dis k, in rough agreement with th e Milky 
Way and FGC 1415 (|Yoa,chim Da,lca,ntonll20f)5^) . How- 
ever, the inclusion of any counter-rotating material in 
the merger could easily break this correlation. For exam- 
ple, the kinematics of FGC 227 indicate that the satel- 
lites which contributed the majority of the baryons to 
the thin disk could not also have deposited the majority 
of the thick disk stars. This particular formation path- 
way allows the scale lengths of the thick and thin disks 
to sometimes decouple, and indeed, the scale lengths of 
FGC 227's thick disk is comparable to, not larger than, 
its thin disk. 

The structural parameters of the thick disks formed in 
the lAbadi ct al. (2003) and Brook et al. (2004, 2005) 
simulations are in moderate agreement with our results. 
However, direct comparisons are difhcult because the 
simulated galaxies tend to be more massive than the 
galaxies in our sample and also host large bulge compo- 
nents. The simulated thick disks do seem to match the 
observed trends of scale height ratios (Figure and lu- 
minosity ratios f Figure [ 191). Howev er, the scale length 
ratio found in the lBrook et alJ l)2005j) simulation is fairly 



low (Figure possibly due to the fact that their thick 
disk stars are formed directly from the gas during merg- 
ers, increasing the likelihood that the thick and thin disk 
stars will share similar scale lengths and kinematics. It 
is also difficult to compare our 2-d decompositions with 
analysis of simulations that can separate stellar popula- 
tions based on kinematics. 

5.3. Further Implications 

Given the excellent fit to the body of data on thick 
disks, we now begin to address other implications of the 
accretion scenario developed above. 

5.3.1. Old Low Mass Galaxies 

In hierarchical galaxy formation models, small scale 
structure collapses first, suggesting that low mass galax- 
ies should be old. This expectation is in direct conflict 
with observations that low mass galaxies almost always 
have blue colors consistent with young stellar popula- 
tions. This difference is one of the most intractable 
failings of the predictions of semi-analytic models (e.g. 
iBel l et al., 2003; van den Bosch 2002). The existence of 
thick disks that dominate the stellar mass of low mass 
galaxies (Figure|23) solves this conundrum. Our observa- 
tions show that low mass galaxies are indeed dominated 
by an old stellar population, but one that is sufficiently 
old, faint, and diffuse that it has no significant impact 
on the observed colors of the young, high-surface bright- 
ness, star-forming thin disk ft |4.4.4|) . We believe that 
semi-analytic models could be brought into alignment 
with the data if they were to include both the locking up 
of material into a diffuse thick disk and the suppression 
of star formation efficiencies in low mass disks due to 
their lying ent irely below the Kennicutt star form ation 
threshold re.g.. lVerde et al.l2002trDalcanton et al.l2 0041. 

5.3.2. Abundance patterns and the timing of thick and 
thin disk formation 

Studies of a-element abundances of the Milky Way 
have suggested that star formation in the thick disk took 
place over several gigayears fe.g.. IBensbv et al . 2004a). 
The abundances in thick disk stars show a flat plateau 
at high [a/Fe] that extends to [Fe/H]'^ — 0.3, indicating 
that thick disk stars enriched quickly to relatively high 
metallicity, before Type la supernovae became preva- 
lent. At larger iron abundances (—0.3 <[Fe/H]< 0), 
however, the a abundance declines linearly, suggesting 
that star formation in the thick disk was sufficiently ex- 
tended (^1 — 3 Gyr) that enrichment from Type la's 
became important. The abundances of thin disk stars 
show similar, parallel behavior, but the plateau does not 
extend to equally high metallicities, indicating that early 
star formation in the thin disk was not nearly as rapid 
as in the thick disk. The abundance patterns of thick 
and thin disk stars therefore follow parallel but distinct 
sequences on the [a/Fe] vs [Fe/H] pl ane, with s ignificant 
overlap in [Fe/H] (rnost recently .Bensby et al. .2004bl: 
IMishenina eiralll200l ICxratton et alJl2003D . 

The above sequence of events is typically taken as evi- 
dence that the thick disk formed from violent heating of 
a previous thinner disk. However, it may be possible to 
accommodate the abundance data in the accretion sce- 
nario as well. First, the rapid enrichment of future thick 
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disk stars can easily occur in the pre-galactic fragments. 
These mini-halos should be dense, leading to high gas 
densities and star formation rates, which would produce 
the necessary fast enrichment. While we have hypoth- 
esized above that supernova-driven winds will truncate 
star formation in the lower mass progenitors, some mini- 
halos will have sufficient mass to retain gas for longer 
periods of time, allowing stars to form over sufficiently 
long timescales to produce both the drop in [a/Fc] and 
the enrichment of some thick disk stars to near solar 
metallicities. 

The expected timescales for this scenario are compat- 
ible with the observational constraints. Theory suggests 
that the epoch of thick disk assembly should correspond 
to the period of rapid mass accretion seen in simulations 
at 2: > 3, or tiookback > 11 Gyr fe.g.. .Zhao et al. 200^. 
Observationally, the relative abundance of [Eu/Ba] indi- 
cates that th ick disk stars were form ed on a timescale 
of 1-1.5 Gyr ijMashonkina et ;ini200.lD . which makes the 
theoretical expectation consistent with the age of the uni- 
verse determined from WMAP. 

In addition to the short star formation timescale for 
thick disk stars, the accretion scenario can produce a 
long timescale for formation of thin disk stars. After the 
pre-galactic fragments merge into a disk, the gas that 
forms the thin disk gradually converts into stars. The 
timescale of this conversion is controlled primarily by 
the gas surface density. In general, this timescale should 
be much longer in the disk than in the pre-galactic frag- 
ments, because the gas is spread over much larger areas, 
leading to lower gas densities and longer star formation 
timescales. The difference in timescales for thick and 
thin disk star formation could lead to the appearance of 
a "delay" between the formation of the two components. 
However, as accretion of both gas and stars would be on- 
going from early times until z ~ 3, some genuinely old 
th ick disk stars woul d be allowed to form (see discussion 
in lAbadi et aDIMiall . 

The accretion scenario also provides a mechanism for 
producing thin disk stars with lower a-abundances than 
thick disks stars at the same mctallicity. Because the thin 
disk assembles from gas that had not been consumed by 
star formation in pre-galactic fragments, it is possible 
for the gas to initially have lower mean metallicity than 
the thick disk stars that were accreted. The gas may 
come from larger radii within individual mini-halos, and 
thus be less enriched. It may also come from fragments 
that have never formed stars, or from cold flow accretion 
directly. Thus, accretion may allow the youngest thin 
disk stars to be sufficiently metal-poor that they overlap 
the metallicities of thick disk stars. 

The one significant trouble spot is the thinness of the 
observed [a/Fe] vs [Fe/H] relation for thick disk stars. 
If the thick disk formed from assembly of many differ- 
ent sub-units of different masses, lifetimes, and gas rich- 
nesses, then one might expect large variations in the de- 
gree of a-enhancement in the accreted stars. On the 
other hand, the potential discrepancy might not be as 
severe as one might initially believe. If supernova-driven 
winds truncate star formation in low mass sub-units, 
then only the most massive mini-halos contribute stars to 
the high metallicity ([Fe/H]> —0.3) thick disk, since they 
are the only precursors that could hold gas long enough 
to allow significant Type la enrichment. Massive halos 



are rarer than low mass halos, and thus a relatively small 
number of halos may dominate the metal rich end of the 
thick disk population, much in the way that galax- 
ies dominate the luminosity density of the local universe. 
These disrupted sa tellites may a lso segregate to different 
radii, as seen in the iAbadi et alJ |2003) simulations, such 
that a sample at the solar circle is dominated by an even 
smaller number of massive satellites. More detailed sim- 
ulations are needed to evaluate the size of this possible 
discrepancy. 

5.3.3. Pre- Enrichment of Thin Disks 

Chemical abundance data on stars within the Milky 
Way has led to the conclus ion that the thin disk may have 
been "pre-enriched" (e.g..lCaimni 3l2OO0tlChia"r)pini et alJ 
IToorHPagel fc Tautvaisienelll995D . Such pre-enrichment 
naturally explains the lack of truly metal poor stars in 
the thin disk as well as the under-abundance of more 
moderately metal-poor stars (i.e. the "G-Dwarf" prob- 
lem). In the scenario we have explored here, the gas from 
which thin disk stars form was originally associated with 
the thick disk, and thus will have been enriched while still 
in pre-galactic fragme nts. While this idea has been sug- 
gested before (e.s.. iBrook et al.l i2005'l. the universality 
of thick disks suggests that it is probably a wide-spread, 
phenomena. 

5.3.4. Producing the Mass- Metallicity Relationship in 

Disks 

Another attractive feature of the satellite accretion 
model is that it facilitates creating the mass-metallicity 
relationship in disks. Many authors have argued that 
the lower metallicities and effective yields seen in low 
mass galaxies is due to the onset of supernova-driven 
winds at the mass scale where the metallicity begins to 
falljK ~ 120 km s ~\ or Mw„„„ < 3 x 10i°Mf^; e.g . 
rxarnettll2002: .TremTmt.i et al]l2004 iDekel fc Wodl2003t 
Kauffmann et al.l20n3|) . However, simulations of gas out- 
flow find that it i s quite difficult to d rive coherent winds 
at these masses l|Ferrara fc ToLstovl [2000') . particularly 
given the low star f ormation rates typical of low mass 
disks fe.g. lHunter fc Elmegreen 2004.') . 

As an alternative, the satellite accretion model sug- 
gests that non-negligible star formation took place in 
lower mass sub-units. These pre-galactic fragments had 
much lower escape velocities, and probably had higher 
gas surface densities due to not yet being organized into 
a coherent rotating disk. Thus, the sub-units are a more 
natural environment for driving winds, given their low 
escape velocities and likely high star formation rates. 
The origin of the observed mass-metallicity relation may 
therefore lie not so much in the disks themselves, but in 
the sub-units from which they assembled. 

6. CONCLUSIONS 

We fit thin and thick disk components to a sample 
of 34 late-type edge-on spiral galaxies. Our thick disk 
components are very similar to previously detected thick 
disk systems and the MW thick disk, suggesting they are 
a remnant stellar population left over from early stages of 
galaxy formation. In lower mass galaxies {Vc < 100 km 
s^^), the thick disk is the dominant component in both 
luminosity and stellar mass. For higher mass galaxies. 
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the thick disk is a minor component, and is analogous to 
the thick disks found in the Milky Way and other higher 
mass galaxies. In particular, we find: 

• Thick disks have a scale height ~ 2 times larger 
than thin disks 

• Thick disks have systematically larger scale lengths 
than thin disks 

• In low mass galaxies, the thick disk can dominate 
the total i?-band luminosity 

• The thick disk comprises 5-40% of the total bary- 
onic mass of our galaxies 

We combine these results with the findings of other 
studies of thick disks to analyze possible thick disk 
formation scenarios. In particular, we include result s 
from thick disk kinematics l)Yoachim fc Dalcantonl2d05() , 
studies of resolved stell ar po pulations in thick disks 
l|Seth et all 120051 iMouldl 120051) and simulations which 
form thick disks f Brook et aLll2004t iBrook et all 120051 
lAbadi e t al. 2003). Overall, we find that models where 
the thick disk forms from a kinematically heated thin 
disk is not supported by the data. Instead, our results 
favor models where thick disk stars formed in galactic 
sub-units before merging to create the final galaxy. 



We consider a hierarchical galaxy formation scenario 
where galaxies form through a series of mergers where 
sub-units deposit both stars and gas. Any stellar com- 
ponent in the sub-units end up in the thick disk, while 
gas cools and forms a thin disk. We find that the low 
mass galaxies in our sample must have formed from sub- 
units that had a higher stellar mass fraction than those 
that formed higher mass galaxies. We can explain this re- 
sult if low mass sub-units (which go on to form low mass 
galaxies) are more susceptible to SN-induced blowout, 
leaving them with a higher stellar to gas mass fraction. A 
mass-dependent blowout scenario is consistent with other 
general observations of disk galaxies, such as the mass- 
metallicity relation and the chemical pre-enrichment of 
the MW thin disk. 
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Single disk fits for the sample galaxies. 
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1 These fits use EQuationsl2landl3lw ith N — 1 (i.e. a sech ^ vertical profile). Peak edge-on surface brightnesses have not been corrected for inclination. 
When available, distances taken from l^arachentsev et al.ll200Cl) . Otherwise, wc have used the recessional velocity corrected for Local Group infall to 
the Virgo cluster (LEDA). Throughout, we assume Hq — 70 km Mpc"i. 



Yoachim & Dalcanton 



TABLE 4 

Vertical light profiles for the two disk models used in fitting R-band structural 

parameters. 



Thin disk model Thick disk model Notes 

sech^(z/zo) sech^(z/2o) convolved w/1" FWHM Gaussian 

sech^(z/2o) sech^(2/2o) midplane masked 

sech^(2;/zo) sech^(z/2o) 

sech(2/zo) sech^(2/zo) 

sech^(2/2o) sech(2/zo) 

sech{2/2o) sech(2/zo) 



Structure of Thick Disks 



TABLE 5 

Vertical light profiles for the two disk models used in fitting i?-BAND structural parameters. 
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